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ABSTRACT 

High resolution spectroscopic observations with FUSE and HST STIS of atomic and molecular velocity strat- 
ification in the nebular outflow of M27 challenge models for the abundance kinematics in planetary nebulae. 
The simple picture of a very high speed (~ 1000 km s ), high ionization, radiation driven stellar wind sur- 
rounded by a slower (^10 km s ) mostly molecular outflow, with low ionization and neutral atomic species 
residing at the wind interaction interface, is not supported by the M27 data. We find no evidence for a high 
speed radiation driven wind. Instead there is a fast (33 - 65 km s ) low ionization zone, surrounding a slower 
(< 33 km s ) high ionization zone and, at the transition velocity (33 km s ), vibrationally excited H2 is 
intermixed with a predominately neutral atomic medium. The ground state H2 ro-vibrational population shows 
detectable absorption from /" < 15 and v" < 3. Far-UV continuum fluorescence of H2 is not detected, but Ly- 
man a (Lya) fluorescence is present. We also find the diffuse nebular medium to be inhospitable to molecules 

and dust. Maintaining the modest equilibrium abundance of H2 ( jmm -C 1) in the diffuse nebular medium 

requires a source of H2, mostly likely the clumpy nebular medium. The stellar spectral energy distribution 
shows no signs of reddening (E(B— V) < 0.01), but paradoxically measurements of Ha/Hp reddening found in 
the literature, and verified here using the APO DIS, indicate E(B— V) ~ 0.1. We argue the apparent enhance- 
ment of Ha/HP in the absence of dust may result from a two step process of H2 ionization by Lyman continuum 
(Lyc) photons followed by dissociative recombination (H2 + y — > H^ + e~ — > H (Is) +H(nl)), which ultimately 
produces fluorescence of Ha and Lya. In the optically thin limit at the inferred radius of the velocity transition 
we find dissociation of H2 by stellar Lyc photons is an order of magnitude more efficient than spontaneous 
dissociation by far-UV photons. We suggest that the importance of this H2 destruction process in H II regions 
has been overlooked. 

Subject headings: atomic processes — ISM: abundances — (ISM:) dust, extinction — (ISM:) planetary nebu- 
lae: general — (ISM:) planetary nebulae: individual (NGC 6853) — line: identification — 
line: profiles — molecular processes — plasmas — (stars:) circumstellar matter — (stars:) 
white dwarfs — ultraviolet: ISM — ultraviolet: stars 



1. INTRODUCTION 

M27 (NGC 6853, the Dumbbell) exhibits a bi-polar mor- 
phology, often associated with molecular hydro gen (H2) in- 
frared emission in planetary nebulae (PNe) dKastner et al.l 
1996). High resolution spectra of the hot central star (CS), 
acquired with the Far Ultraviolet Spectroscopic Explorer 
(FUSE), have revealed an unusually rich set of narrow H2 
absorption features spanning the entire spectral bandpass, an 
indication that the molecule is vibrationally excited. FUSE 
carries no on-board source for wavelength calibrations and 
consequently M27 has been observed numerous times for this 
purpose, resulting in a high signal-t o-noise data set (see Fig- 
ureQ]and McCandlis s & Krukll2007l) . The usefulness of these 
lines for wavelength calibration aside, we are presented with 
an interesting puzzle. What physical processes excite H2 in 
PNe and how does it survive in this high temperature and 
highly ionized environment? 

At first glance it seems surprising to find H2 in PNe at all. 
For a typical electron density of 300 cm~ 3 and temperature 
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of 10,000 K, we estimate an e-fold lifetime w 230 years, as- 
su ming electron imp act dissociation of H2 at the rate given 
bv lMartin etakl d!998l y e - = 4.6 x 10~ 13 cm 3 s" 1 ). Nev- 
ertheless there is a whole class of PNe with bi-polar mor- 
phology i n which infrared H2 emission is the defining char - 
acteristic dKastner et all 1 1 996t IZuckerman & Gatlevl Tl988b . 
Zuckerman & Gatleyj concluded that the infrared H2 emission 
in M27 is shock excited, based on IR spectroscopic diagnos- 
tics. This result is considered somewhat surprising because 
the close presence of a hot star suggests that far-UV contin- 
uum ex cited fluorescence of t he H2 may be important. For in- 
stance, iNatta & Hollenbachl (Q998) calculate that continuum 
pumped fluorescence processes dominate thermal processes 
in PNe evolutionary models at la te ti mes > 5,000 yrs. 

iHerald & Bianchil(l2002ll2004 and lDinerstein etaTJ(l2004 
have reported excited H2 absorption features in FUSE spectra 
of a handful of CSPN located in the Gal axy and Large Magel- 
lanic Cloud. iDinerstein et all (12004 and lSterling et all (l2005h 
have used FUSE to examine PNe with strong, extended H2 in- 
frared (IR) emission. They find the detection of IR emission 
is no guarantee for finding excited H2 in absorption and con- 
clude the molecular material in these systems is clumped. 

Many authors have discussed the clumped structures and 
dense knots observed in PNe, along wit h evidence for their 



origin and chemical composition (c.f. lO'Dell et all I2002L 
2003 1 : iHuggins et all l2002t iBachiller et al.l 120001: ICox et all 
1998; Meaburn & LopeJ l 19931: iReav & Athertonlll985L 
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FIG. 1 . — Flux of the orbital coadds for the eight spectral segments taken through the HIRS spectrograph aperture. Individual segments are labeled. Overplotted 
is the model SED without extinction (upper line) and with E(B— V) = 0.005 mag (slightly lower line) and E(B—V) = 0.05 (lowest line) for a standard galactic 
extinction curve. E(B—V) = 0.005 mag provides the best match to the SED in the SiC channels. 



references therein). These structures are isolated photodisso- 
ciation regions (PDR) immersed in H II regions, and it is nat- 
ural to expect them to be reservoirs of molecular material in 
PNe, regardless of the mechanism that causes them to form. 
Capriotti (1973) studied the dynamical evolution of a radia- 
tion bounded PN and found that dense neutral globules form 
as the result of a gradually weaking radiation field in an out- 
ward propagating ionization front. A Rayleigh-Taylor like in- 
stability develops and eventually an optically thick spike of 
neutral material becomes separated from the front and forms 
a high density globule where molecules can presumably form. 
The formation of H2 at later stages, on dust grains and in 
the high electron density environment via the H + H~ — > H2 
+ e~ react ion, has been modeled by lAleman & Gruenwaldl 
(12001 and lNatta & Hollenbachl (119981) . IWilriamsId 19991) has 
di scussed a shadow ing instability, somewhat similar to that 
of ICapriottH (119731) . where small inhomogeneities perturb a 
passing supersonic ionization front, causing a corrugation that 
produces a large ne utral density contrast further downstream. 
iDvson et al.1 d!989l) suggested molecular knots in PN are relic 
SiO maser spots, which originate in the atmosphere of the 
asymptotic giant branch (AGB) progenitor. Others have pro- 
posed that molecules may also originate from relic planetary 
material, either accreted and then ejected or swept up during 
the AGB phase preceding the formation of the nebula (c.f. 
Wesson & Liu! 120041 iRvbicki & Denisl 1200 It ISiess & Liviol 
19991; iLivio & Soker|[l983l) . ISokerl d!999l) has discussed the 
si gnatures of surviving Uranus/Neptune-like planets in PNe. 

iRedman et al.1 (120031) have emphasized the importance of 
molecular observations as a means to establish the evolution- 



ary history of the clumps, which should depend on when and 
where they formed during the AGB — > PN transition. The 
expectation is for dense clumps formed within an AGB at- 
mosphere to have more complex molecules, due to higher 
dust extinction and molecular shielding, than clumps formed 
later in the PN phase when the stellar radiation field is harder 
and the overall density lower. The abundance of atomic and 
molecular species with respect to H2 is of fundamental impor- 
tance to assessing the predictions for the chemical evolution 
in PNe clumps. Observations of atomic and molecular excita- 
tion stratification in the outflow provide a fossil record of the 
AGB — » PN transition containing clues to the origins of the 
clumps. 

Towards this end we explore high resolution far-ultraviolet 
absorption spectroscopy, provided by FUSE and the Space 
Telescope Imaging Spectrograph (STIS) on board the Hubble 
Space Telescope {HST), to reveal the nebular outflow kine- 
matics imprinted on the line profiles of a wide variety of 
molecular and atomic species, including H2, CO, H I, C I - 
IV, N I - III, O I, O VI, Si II - IV, P II - V, S II - IV, Ar I - II and 
Fe II - III. Absorption spectroscopy samples material directly 
along the line-of-sight. We will sometimes refer to it as the 
diffuse nebular medium to distinguish it from the extended 
clumpy medium, offset from the direct line-of-sight, which 
undoubtedly differs in physical and chemical composition. 

Our main objective is to quantify the H2 excitation state, its 
abundance with respect to H I, and its outflow velocity with 
respect to other emitting and absorbing atomic and molecu- 
lar species. We use this information, along with the spectral 
energy distribution (SED) of the central star, to constrain H2 
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FIG. 2.— Normalized intensity profiles of H a , [O III] ^.5007 and [N II] 
X 6584 as a fun ction of velocity. H a and [N II] X 65 84 were digitized from 
IMeaburrj J2005l) and [O III] ^.5007 was digitized from lMeaburn et all (199% . 

formation and destruction processes in the nebula. 

The analysis is supported by a variety of ancillary data 
on M27. Longslit far-ultraviolet spectroscopy from a 
JHU/NASA sounding rocket (36.136 UG) provides an upper 
limit to the H2 continuum fluorescence. Longslit optical spec- 
trophotometry with the Double Imaging Spectrograph (DIS) 
at Apache Point Observatory (APO) yields the absolute flux 
of the CS and the surrounding nebular Balmer line emission. 
Dwingeloo Survey 21 cm data helps constrain the velocity 
structure of H I. We begin by reviewing previous observations 
and physical properties of the nebula followed by a discus- 
sion of the various data sets. We then present the analysis, 
discussion, conclusions and suggestions for future work. 

2. PHYSICAL PROPERTIES AND PREVIOUS OBSERVATIONS OF 

M27 

M27 h as an elliptical sh ape w 8' x 5' as seen in visual pho- 
tographs (lBurnhamlll978l) . Deep narrow band CC D imaging 
of Ha + [N II] and [O III] bv lPapamastorakis et alj d 1993b has 
shown a faint halo of these emissions extending to 17'. Its 
bi-polar morphology is manifest in Ha images as t he "D umb- 
bell" and in the 2.12 /jm image of lKastner etall 0996) as a 
clumpy "bow-tie" aligned with the semi-minor axis. The ends 
of the bow-tie are ~ 3' wide, the waist ~ 1' and the len gth w 
6' (c.f.lKastner et alJ|1996UZuckerman & Gafleyll 19881) . 

Napiwotzki (1999) determined the DAO CS parameters 
from NLTE model atmosphere analysis and evolutionary con- 
siderations (T = 108,600 ± 6800 K, log(g) = 6.7 ± 0.23, 
log(He/H) by number of -1.12 dex, and M* = 0.56 ± 0.01 
M ). Astrometric observations of the CS by iBenedict et alj 
(2003) produced a distance of d = 417+gj pc. At this distance 
the w 8'x 5' ellipse is 1 pc x 0.6 pc an d the 1025" diame- 
ter of the halo is 2.1 pc. Benedic t et alj give other physical 
parameters as well, V =13.98±0.03, a total extinction Ay = 



0.3±0.06 and a stellar radius R t = 0.055±0.02 R . They find 
a bolometric magnitude M* bol = -1.67±0.37, which yields a 
luminosity for the central star of L t = 366 L@. The radial 
velocity of t he central star in the heliocentric system was de- 
te rmined by I Wilson! (119531) to be V sys = -42 ± 6 km s -1 . 

Barker ( 1984) has reported on the analysis of line emission 
from the nebula as derived from both ground based and In- 
ternational Ultraviolet Explorer (IUE) spectra. The emission 
is typical of an H II region with an electron temperature of 
« 10,000 K, an electron density of w 300 ± 100 cm" 3 , and 
an elevated metallicity with respect to sola r in CNO. These 
param eters agree well with those found by lHawley & Millet 
i[T978|) in a similar ground based study of the nebula. iBarkerl 
( 1984) adopts a line reddening parameter of c = 0.17 (c is de- 



fined in § 14.2. 11 1. The literature reveals a r ange of nebular line 
reddening parame t ers, 0.03< c < 0.18 (c.f.lMillerfl973tlKaleri 
[l976t ICahnlll97l lBarkerffi984l ~ ICiardullo et al.lll999l) . This 
translates to 0.02< E (B-V) < 0.12, (using clE(B-V ) = 1.5, 
ICiardullo et aTTl 999. Figure 4). iPottasch et all (1 19771) derived 
E{B— V) = 0.10 ± 0.04 by removing a small (some would 
say imperceptible) 2175 A bump from an IUE spectrum. This 
valu e is nearly the s ame a s that given by lHamset ail (119971) 
and IBenedict et al. (2003), based on trigonometric distance 
determinations an d absolute visual magnitude considerations. 
ICahn et alJ (1 19921) report c r = 0.04 derived from the ratio of 5 
GHz radio flux density to HP absolute flux. 

Velocity and positio n resolved emission line spectroscopy 
bv lGoudis etalJ(ll978l) showed faint [O I] and [N II] profiles 
(beam size of 83") taken near the CS to be double peaked 
symmetrically around the radial velocity of the system. In 
the outer regions the profiles converged to single peaks at the 
systemic velocity. The bright [O III] lines (beam size of 30") 
showed similar behavior but with a lower double peaked split- 
ting. These observations are consistent with the presence of 
two nested shells expanding w ith projected v e locitie s of « 
33 and 15 km s respectively. iMeaburn et al.1 (1 19921) used a 
fiber optic image dissector fed into an echelle spectrograph to 
closely examine the O III velocity structure in the immediate 
vicinity of the CS (« 36" North-South), and found quadruple 
peaked line profiles consistent with gas expanding a t veloci- 
ties of 31 km s and 12 km s -1 . IMeaburrj (120051) showed 
that Ha exhibits a filled-in asymmetric profile with a half- 
width half-maxium of w 37 km s , in contrast to the double 
peaked emission exhibited by [N II] A.6584 in the vicinity of 
the CS. He II A.6560 shows a much narrower profile and is 
consistent with turbulent broadening alone. Example s of the 
[N II] A.6584 and Ha profiles from IMeaburn! (120051 Figure 
3) are reproduced in the top and bottom panels of Figure [2] 
T he middle panel of Fig ure [2] shows O III A.5007 taken from 
of IMeaburn et al.1 (1 19921 Figure 6d). Note the base of these 
emission features extend to ± 60 km s~ 1 . 

Observations bv lHuggins etal.1 {l996) of the CO(2-l) 230 
GHz emission lines at a position 68"W and 63"S from the 
CS also exhibit similar double peaked line profiles, with a to- 
tal separation w 30 km s _1 . At this position they estimate 
a CO column density of 6.9 x 10 15 cm~ 2 (assuming LTE 
with 5 < r CT <150K) . At a position 10" due West of the CS, 
iBachiller etaT] (2000) show lines split by 53 ± 1 km s _1 (their 
Figure 4), yielding an expansion velocity of w 27 km s _I 
with symmetry about a heliocentric velocity of ^43 km s _1 (at 
M27, Vh e i = Vi sr - 17.8 km s~'), in good agreement with the 
Wilson systemic velocity. They present a map of the region 
(beam size 12"), showing the CO emission to be clumpy and 
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FIG. 3. — Top - raw longslit far-UV spectrum from sounding rocket experiment 36.136 UG. Bottom - total nebular brightness per A summed over the whole 
slit, excluding the continuum source. The red line is the estimated detector background. There is no significant dust scattering or H2 continuum fluoresence 
emission above this background. The rise longward of Lya is likely caused by 2s - Is, 2-photon emission of recombining hydrogen. The dotted blue line shows 
the expected brightness for a column of N(2s) = 4 x 10 8 cm~ 2 . 



more-or-less coinciden t with the mol e cular h y drogen infrared 
emissi on observed by iKastner et alj 0996). Bachill er et alj 
(2000) present a picture of an ionized central region with an 
electron density of about n e ~ 100 cm~ 3 surrounded by a 
ring of molecular clumps with densities ~ 10 4 cm~ 3 under- 
going photodissociation. They suggest the clumps are sim- 
ilar to the "cometary" sha ped features observed in great de- 
tail in the Helix nebula (c.f . lO'Dell et al.ll2002t iHuggins et ail 
l2002t iMeixner et al.ll2005l) . although the cometary morphol- 
ogy is not as prominent in th e Dumbbell as it is in the He- 
lix. iMeaburn & Lopej 993) first noted (see also HST im- 
ages bv lO'Dell et alj|2002l) that some clumps appear as dark 
knots against the very bright O III emission that surrounds the 
c entral region. 

iMcCandlissl (2001) gave a preliminary analysis of molecu- 
lar hydrogen absorption in the FUSE spectra, identifing two 
absorption systems, one blueward of the CS radial velocity 
and one redward. The blueward component has lines originat- 
ing from an excited electronic ground state (X with ro- 
tational levels (J") up to 15 and vibrational levels (v") at least 
as high as 3. Many lines originating from v" = 1 and 2 are lo- 
cated longward of the ground state bandhead (V - v" = — 0) 
at 1 108.12 A, and as such are unambiguous markers of "hot" 
H2. In contrast, the redward component is much "cooler," 
showing no unusual excitation (J" mx « 4, v" = 0). It is associ- 
ated with non-nebular foreground gas. Here we will refine the 
original analysis, which was based on LWRS data acquired 
under FUSE PI Team Program PI 04. 

Lastly we note, Lya fluorescence of H ? has recent l y bee n 
found in the PNe M27 and NGC 3132 bv lLupu et al] (I2006T) . 
The fluorescence is a direct consequence of excited molecular 
hydrogen, as the mechanism requires a significant population 



in the v" = 2 levels in the presence of strong Lya emission. 
Herald reports (private communication) that FUSE observa- 
tions of the CS of NGC 3132 also exhibit excited H2 absorp- 
tion lines with a roughly thermal ro-vibrational distribution of 
wl750K. 

3. DATASETS 
3.1. FUSE 

A description of the data set and processing procedures for 
the FUSE spectra of the CS (M AST object ID GCRV12336 ) 
used in this study can be found in lMcCandhss & Krukl(l2007l) . 
For a description of the detectors, channel alignment issues, 
systems nom enclature an d othe r asp ects of the FUSE in - 
strument, see iMoos et al.l (l2000h and ISahnow et all (120001) . 
Briefly, the CS was observed numerous times for observatory 
wavelength calibration purposes. High signal-to-noise spec- 
tra, acquired in time tag mode through the high resolution slits 
(HIRS) of the eight different channel segments have been as- 
sembled. Data for each extracted channel segment consists 
of three 1-D arrays: wavelength (A), flux (ergs cm~ 2 s _1 
A -1 ), and an estimated statistical error (in flux units). Fig- 
ure Q] shows the eight individual channel segments and their 
spectral range for all the coadded HIRS spectra. The high 
density of the H2 features is apparent. 

In the low sensitivity SiC channels the continuum signal- 
to-noise is between 10-25 per pixel, while in the higher 
sensitivity LiF channels it is 20-55. These signal-to-noise 
ratios are purely statistical and do not account for systemic 
errors, such as detector fixed pattern noise. The resolution of 
the spectra changes slightly as a function of wavelengt h fo r 
each channel segment. In modeling the H2 absorption (§ 14. U . 
we find that a gaussian convolution kernel with a full width at 
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FIG. 4. — Log-log plot of the M27 c entral star SEP from 900 - 20000 A, as measured by FUSE, IUE, APO and the optical photometry of Tvle nda et all 
1 1991) (o), i Ciardullo et al. 1999) (□) and Benedict et al. 12003) (+). Overplotted in black, purple and red are the respective unextinguished, lightly extinguished 
{E(B— V) = 0.005) and heavily extinguished (E(B— V) = 0.05) stellar continuum models. Break between 2000 and 3000 A indicates a gap in the stellar model. 



half maximum of 0.056 A at 1000 A (spectral resolution R w 
18,000, velocity resolution 8V « 17 km s _1 ) provides a good 
match to the unresolved absorption features throughout most 
of the FUSE bandpass. 

Close comparison of the wavelength registration of overlap- 
ping segments reveals isolated regions, a few A in length, of 
slight spectral mismatch (~ a fraction of a resolution element) 
in the wavelength solutions. Consequently, combining all the 
spectra into one master spectrum will result in a loss of reso- 
lution. However, treating each channel/segment individually 
increases the bookkeeping associated with the data analysis. 
Further, because fixed pattern noise tends to dominate when 
the signal-to-noise is high, there is little additional informa- 
tion to be gained in analyzing a low signal-to-noise data set 
when high signal-to-noise is available. For these reasons we 
elected to form two spectra, each of which covers the 900 - 
1 190 bandpass contiguously, using the following procedure. 

The flux and error arrays for the LiFla and LiFlb seg- 
ments were interpolated onto a common linear wavelength 
scale with a 0.013 A bin, covering 900 - 1 190 A. The empty 
wavelength regions, being most of the short wavelength re- 
gion from 1000 A down to the 900 A and the short gap re- 
gion in between LiFla and LiFlb, were filled in with most 
of SiC2a and a small portion of SiC2b respectively. We refer 
to this spectrum as sl2. The LiF2b, LiF2a, SiCla and SiClb 
segments were merged similarily into a spectrum, s21. Ab- 
sorption line analyses were carried out using the composite 
sl2 and s21 spectra (§ 14. U . 5 



5 The s!2 and s21 spectra are available through the H200IS website 
iMcCandliss 2003) along with the original complete Ml 07 data product 
(LWRS, MDRS and HIRS) processed by Jeffrey Kruk (private communica- 
tion). 



3.2. Far-UV Longslit Observations of M27 from JHU/NASA 
Sounding Rocket 36.136 UG 

JHU/NASA sounding rocket mission 36.136 UG was 
launched from White Sands Missile Range Launch Complex- 
36 on 14 June 1999 at 01:40 MST Its purpose was to ob- 
serve the hot CS of M27, provide spatial information on the 
excitation state of H2 in the nebula and investigate its dust 
scattering properties. The science instrument was a 40 cm di- 
ameter Dall-Kirkham telescope with SiC over-coated Al mir- 
rors, feeding a Rowland circle spectrograph with a 200" x 12" 
longslit, a 900 - 1400 A bandpass and an inverse linear dis- 
persion of ~ 20 A mm' 1 . The basic config urati on has been 
describ ed by IMcCandliss et alJ (11994 120001) and iBurgh et all 
(|2001l) and two sim i lar missions using th is payload have flown 
(Burg h et alJ l2002t iFrance et all 120041) . This was the first 
flight for a newly reconfigured sp ectrograph with a holo graph- 
ically corrected concave grating (IMcCandliss et al.ll2001l) . for 
improving the spatial resolution (3" spacecraft pointing lim- 
ited) and a high QE KBr coated micro-channel plate with 
double-delay line anode (~ 25 ^m reso lution element), simi- 
lar to that used by iMcPhate et all (119991) . 

The target was acquired at w T+150 s, at which time point- 
ing control was passed to a ground based operator, who used 
a live video downlink of the nebular field to make real time 
pointing maneuvers. At w T+500 s the telescope was sealed 
prior to reentering the atmosphere. The payload was recov- 
ered and post flight calibrations were made to secure an abso- 
lute calibration. 

The pointing corrected two dimensional spectrum is shown 
in the top panel of Figure|3] C II A.1335, C III A.977 and XI 175, 
He II+N II A.1085, N III A.A.989 - 991 and H I Lyman emission 
lines are evident along with the continuum spectrum of the 
central star. The bottom panel shows the extracted nebular 
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FIG . 5 . — Long slit profiles of Hcc/HP, HpVHy, and Hp/H8. Intrinsic line ratio for a gas with a temperature of 1 0000 K and a density of 1 00 cm -3 is overplotted 
in black. The data have been binned up by two 0"4 pixels. Average values in the SW and NE region, delimited with the vertical dashed lines, are overplotted in 
red. Error bars are displayed with respect to the average in the SW region. All the ratios are consistent with zero extinction except for the NE Ha/Hp ratio. This 
result suggests that a mechanism, other than extinction by dust, is causing the Ha/Hp ratio to deviate from the expectations of recombination. 



spectrum integrated over the length of the slit, excluding the 
stellar continuum. The strong Lyoc emission feature that fills 
the slit is dominated by geocoronal emission. There is a hint 
of H I 2s - Is two photon emission starting at Lyoc and ris- 
ing towards 1400 A. The dotted blue line shows the expected 
brightness for an column of N (H(2s)) = 4 x 10 8 cm -2 (e.g. 
Nussbaumer & Schmutz 1984). There is no evidence for con- 
tinuum pumped fluorescence of H2 or dust scattered stellar 
continuum down to the detector background limit of » 5 x 
A -1 arcsec -2 indicated by the red line. 
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3.3. APO DIS Longslit Observations ofBalmer Line Ratios 

in M27 



It is curious that lHawley & Milier! (119781) and iBarkeri 
d 19841) both cited problems with their Ha /H$ ra ti os be - 
ing higher than the 2.9 ratio determined by iMillerl d!973l) . 
lHawley & Milled (1 19781) suggested some unidentified system- 
atic er ror was affecting the ratio at the 10% level. Barker 
( 1984) reached a similar conclusion and went so far as to de- 
crease the Ha/H$ by 13% before applying an extinction cor- 
rection of c = 0.17. Both studies find « 20% variations in 
HaJH$ at different locations within the nebula. 

To more thoroughly investigate this phenomenon we re- 
cently acquired an extensive series of longslit spectral scans of 
the nebula with the Double Imaging Spectrograph (DIS ver- 
sion III) at the ARC 3.5 m telescope at Apache Point Observa- 
tory (APO), during the nights of 1 - 2 July 2006 commencing 
at 01:15 MDT DIS has blue and red spectra channels with 
back-illuminated, 13.5 pixel Marconi CCDs in a 1028 x 
2048 pixel 2 format for recording data. The medium resolu- 



tion gratings have inverse linear dispersions of 1.85 and 2.26 
A pixel -1 respectively. The useful spectral ranges are 3700 - 
5400 A for the blue and 5300 - 9700 A for the red. The plate 
scales are 0"42 pixel -1 for the blue and 0"40 pixel -1 for the 
red. 

The hot sub-dwarf star BP +28 4211 was used as a spec- 
trophotometric standard (Bohl in et all 12001). This star was 
observed through the large 5" x 300" DIS slit for 60 seconds 
at an airmass of 1.18. The M27 CS was observed shortly 
thereafter for 200 seconds through the same aperture at an 
airmass of 1.02. The CS and surrounding nebula were then 
observed through the narrowest 0"9 x 300" DIS slit for 200 
seconds, with the slit held fixed on the CS. The narrowest 
slit provides clean separation of the [N II] A.A.6548, 6584 lines 
from Ha. The position angle of the slit was 35°. Biases, flat 
fields and emission line spectra were recorded at twilight, and 
the data were reduced using custom IDL code. The CS abso- 
lute flux is shown in Figure |4] and will be discussed in § 14.21 
The Balmer line longslit profiles centered on the CS, and ra- 
tioed to //p, are shown in Figure [5] and will be discussed in 
§ 14.2.11 Analysis of the full dataset from the APO run will be 
presented elsewhere. 

3.4. Dwingeloo H 1 21-cm Observation ofM27 

The shapes of atomic hydrogen absorption profiles are com- 
plex, as they result from ensembles of intervening ISM ab- 
sorption systems, separated in velocity space along the line- 
of-sight. Determining the column density and doppler param- 
eter for individual velocity components within these broad 
damped and saturated absorption profiles is difficult without 
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FIG. 6. — Nearest neighbor H I 21 cm emission spectra from the Dwingloo spectral survey atlas. These spectra were used to form a minimum spectrum and a 
distance-weighted mean spectrum. 



a priori information of the line-of-slight velocity structures. 
One way to gain this information is to examine H I 21 cm 
emission data such as found in the Atl as of Galactic Neu- 
tral Hydrogen (lHartmann & Burtonlll997L also known as the 
Dwingeloo survey). The high spectral resolution of the data 
is excellent for locating velocity components, provides an up- 
per limit to the amount of neutral hydrogen in the foreground, 
and serves as a starting point for absorption line modeling. 
The atlas, in galactic coordinates (l,b) with (0.°5) 2 cells, has 
a velocity spacing of 1.03 km s _1 . The positions cover most 
of the sky and the velocity coverage spans -4-50 < v/ sr < 400 
km s _1 in local-standard-of -rest coordinates. The intensity is 
given in antenna temperature per velocity (K (km s ) ), but 
can be converted to a H I column density, under the assump- 
tion of optically thin emission, b y multiplying with a conver - 
sion factor (0.1 82 x 10~ 19 cm' 2 . lHartmann & Burtonll 19971) . 

The galactic coordinates of M27 are / = 60P84, b = -3.°70. 
The neutral portion of the n ebula should be at least as l arge 
as the Ha halo detected by Papamasto rakis etaP (1 19931) (> 
0.°28) and since the Dwingeloo half power beam width is 
slightly larger than the atlas grid, it is reasonable to expect 
some signal from the nebula to appear in the four nearest 
neighbor points, i.e. atlas coordinates (/ = 61?, b = ~3.°5) (/ 
= 60?5, b = -3.°5) (I = 61.°, b = -4.°) and (/ = 60.°5, b = -A?). 
The column density profiles of H I as a function of velocity 
for these four grid points, in heliocentric velocity coordinates, 
are shown in Figure [6] 

To reduce the non-local contributions we constructed a 
nearest neighbor minimum spectrum under the assumption 
that non-local contributions are widespread and uniform. This 
assumption has been checked by examining an area out to 
±1.°5, which showed similar structures, typified by strong 
peaks near km s and smaller peaks near -85 km s . 
The nearest neighbor minimum spectrum was smoothed with 



a five bin boxcar average and subtracted from a similarily 
smoothed spectrum formed from a distance weighted aver- 
age of the same 4 nearest neighbors. This process is shown 
in Figure 13 The resulting subtraction was fit with nine Gaus- 
sian profiles spread between -99 < Vf, e i < 24 km s . The 
integrated emission column densities, rms velocity widths (b 
values), and heliocentric velocities are given in columns 2, 3, 
and 4 of Table |2 

Column 5 shows the column densities that we adopted for 
the model of the Lyman series absorptions as discussed in 
§ 14.31 Absorption line components to the blue of the systemic 
velocity at -42 km s occur within the expansion velocity 
inferred from the optical emission lines and can plausibly be 
associated with the nebular expansion. If the absorption com- 
ponents to the red of -42 km s _1 are associated with the neb- 
ula they would have to be infalling. They are more likely to 
be located in the foreground. 

3.5. STIS 

The E140M spectrum of M27 (o64d07020_xld.fits), taken 
from the Multimission Archive at Space Telescope (MAST), 
was acquired for HST Proposal 8638 (Klaus Werner - PI). 
These data are a high level product consisting of the ex- 
tracted one-dimensional arrays of flux, flux error and wave- 
length for the individual echelle orders. The data were ac- 
quired through the 0"2 x 0."2 spectroscopic aperture for an ex- 
posure time of 2906 s and were reduced with CALSTIS ver- 
sion 2.18. The spectral resolution of the E140M is given in 
the STIS data handbook (version 7.0 lOuijanoll2003l) asfi« 
48,000 (8V = 6.25 km s" 1 ). The intrinsic line profile for the 
0"2x0"2 aperture has a gaussian core with this width, but 
there is non-negligible power in the wings. Comparison of 
the FUSE and STIS wavelength scales for the CS revealed a 
discrepant offset. The establishment of a self consistent wave- 
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FIG. 7. — Top - Minimum spectrum is subtracted from the angular distance weighted mean spectrum to remove wide angle background emission. Vertical axis 
is column density per km s ■ Horizontal axis is heliocentric velocity. Bottom - Resulting spectrum contains both foreground and some background emission. 



Nine component gaussian fit is shown along with the fit residuals (offset by -5 
parameters derived from gaussian fits. 

length scale is discussed in § 14.51 

4. DATA ANALYSIS 

Our primary purpose is to quantify the excitation state of 
the H2 in the diffuse nebular medium. The column densities 
of the individual ro-vibrational states were determined using 
curve-of-growth analysis. We detect essentially imperceptible 
reddening of the stellar SED over a decade of wavelength cov- 
erage, in contrast to the reddening determinations using the 
B aimer decrement method as reported in the literature (§|2j- 
We show new high precision longslit Balmer decrement ra- 
tios (Ha/HP, Hy/Hp\ and H8/Hp) where Hoc/H(3 appears to 
be reddened in certain regions of the nebula, but the higher 
order decrements appear unreddened everywhere. A potential 
resolution of the reddening conundrum will be discussed in 

§m 

We also produce an estimate of the total atomic hydrogen 
column density in the nebula to enable a discussion of H2 for- 
mation and destruction processes. An upper limit is placed on 
the column density of CO in the diffuse medium and an anal- 
ysis of the excitation of C I fine structure lines will provide 
a diagnostic of nebular pressure and insight into the H2 exi- 
tation mechanism. Absorption line profiles of high and low 
ionization species compared to the neutral profiles reveal the 
ionization stratification of the outflow and provide a velocity 
description of the nebular kinematics. 

4.1. Curve-of -Growth for the Molecular Hydrogen 
Absorption Components 

It seems with a first glance at Figure Q] that determining 
the column densities for the ro-vibrational transitions and 
identifying velocity components will be daunting. However, 
the regular spacing of H2 lines arising from common rota- 



10 18 cm 2 ). See Tableffjfor total columns, velocity centroids, and doppler 

tional states, along with the monotonic variation in transition 
strength with increasing upper vibrational level, makes identi- 
fication of discrete velocity components straightforward. The 
high line density becomes an advantage, virtually guarantee- 
ing a few lines for a given rotational state will be unblended. 
This allows column densities to be determined by a straight 
forward measurement of the equivalent widths for use in a 
curve-of-growth analysis. 

We have identified two distinct components with heliocen- 
tric velocities of -75 and -28 ± 2 km s . The blueward 
velocity component is the more highly excited, allowing the 
construction of curves-of-growth for all the ro-vibrational lev- 
els with J" < 11 and v"< 1 of the ground electronic state 
X 6 The red ward component has no extraordinary exci- 
tation with measurable absorption only found for /" < 3 in 
v" = 0. Figure [8] shows a small region around the (v'—v") = 
(4-0) band near 1050A where the "blue" and "red" velocity 
components are identified. Blueward absorptions for (v'-l), 
(v'-2) and (v'-3) in this region are also displayed. 

A semi-autonomous method was developed to measure the 
equivalent widths of the line s, based in part on th e completely 
autonomous method used bv lMcCan dliss (2001). All the lines 
for a given velocity offset and rotational state were located 
within a summation interval initially 13-15 pixels wide. The 
pixels adjacent to either side of the summation interval were 
used to define continuum points, through which a straight line 
was fit to serve as the model for the continuum in the summa- 
tion region. The degree of blending was assessed interactively 

6 For rotational and vibrational quantum numbers in a transition leading to 
an upper electronic state (e.g. either B i Yj l or C'n„) from the ground state 
(X 'Eg ) the convention is to designate the upper state with a single prime (/', 
v') and the lower state with a double prime (/", v"). 




FIG. 8. — Region near the H2 (4-0) band identifying the high excitation velocity component at -33 km s (blue) and the low excitation component at +14 
km s (red) as measured in the CS restframe. FUSE si 2 and s21 spectra are overplotted in orange and green respectively. 



by plotting the spectrum surrounding each individual line in 
a (-150, 50) km s _1 interval and overplotting all the known 
atomic and molecular features. Lines with evidence of blend- 
ing were rejected from the curve-of-growth. Examination of 
the lines for blends provided an opportunity to fine-tune the 
summation interval and continuum placement. The equiva- 
lent width measurements were carried out on both the sl2 
and s21 spectra (see § 13. U . using the same summation and 
continuum intervals. Allowance was made for the slight local 
mismatches in velocity scale that exist between the sl2 and 
s21 spectra, by locating the minimum within the integration 
region and adjusting the interval accordingly. 

Initial curves-of-growth were constructed from the equiva- 
lent width meas urements using the w avelengths and oscillator 
strengths from Ab g rail et alJ (Il993allbl) . We required at least 2 
unblended lines for these curves. Lines from transitions with 
J" > 11 and v" > 2 are detected, but they are weak and too 
few to construct a reliable curve-of-growth. Independent % 2 
fittings were performed for each curve by varying the doppler 
parameter over 2 - lOkms -1 and the logarithm of the column 
density (in cm 2 ) over 13-18 dex. The doppler parameters 
thus derived varied from 4 to 8 km s _1 and the logarithmic 
column densities ranged over w 13.5 - 17 dex. 

Comparing results from the sl2 and s21 measurements re- 
vealed a few disagreements for the doppler parameter (and 
thus column density) derived for the same rotational state. 
This produced an inconsistent modulation in the ortho-para 
ratio, expected to be w 3:1 for the odd to even rotational 
states. The inconsistencies were traced to curves-of-growth 
with the most saturated lines and lowest number of points. 
Restricting the allowed doppler parameters in the fit to either 
6 or 7 km s (the most common values) resolved the col- 
umn density discrepancy and resulted in consistent ortho-para 
modulation. It also reduced the scatter when all the curves- 



TABLE 1 

Derived Nebular H 2 X '£+ Column 
Densities* 
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log/Y(v" = 0) 
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Systematic error = ± 0.1 dex. 
The doppler parameter is 6.5 ± 0.5 km s 



of-growth for the individual rotational states from the sl2 and 
s21 spectral measurements were combined (top panel of Fig- 
ure!*)]). The column densities for the individual ro-vibrational 
levels are listed in Table Q] along with the number of lines 
used for each curve-of-growth. The errors are dominated by 
placement of the continuum. Experiments were performed to 
gauge the effect of systematic continuum offsets on the de- 
rived column densities. We found the column densities were 
consistent to within ± 0.1 in the dex for offsets within the 
local continuum signal-to-noise. 

The total column density for the hot H2 component at -75 
km s _1 is N(H2)=7.9 ± 2 x 10 16 cm 2 with a doppler velocity 
of 6.5 ± 0.5 km s -1 . Curve-of-growth analysis on the cold 
non-nebular component, for which we wish only to identify 
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FIG. 9. — Upper - curves-of-growth for doppler parameters of 6 and 7 km s _I along with the data. A doppler parameter of 6.5 km s _1 was adopted. Lower 
- column densities, normalized by level statistical weight, as a function of the level energy. The v" = 0, < J" <11 columns are shown with + marks, and the 
v" = 1, < J" <11 columns are shown with * marks. Overplotted (solid line) is a single temperature Boltzmann distribution for T =2040 K. There are statically 
significant deviations from the single temperature population, although the higher 7 < J" < 11 states are in agreement. The lower < J" < 6 states appear to 
have a flatter slope indicative of a temperature a factor of two higher as the dashed lines indicate. 

the strength of its aborption features, yielded a total column 



density of 1.3 ± 5 x 10 17 cm 2 with a doppler velocity of 5 
km s~ 1 . The rotational distribution is well approximated with 
a temperature of 200 K. 

In the bottom panel of Figure [9] we show the population 
density of the ro-vibrational levels as a function of level en- 
ergy. We use (+) marks to indicate the columns for 7" in the 
v" = state and (*) marks to indicate the columns for J" in 
the v" = 1 state. The straight solid line is the best fit single 
temperature model (T = 2040 K), assuming a Boltzmann dis- 
tribution of level populations. This model is not successful, 
as there are statistically significant deviations from a single 
temperature Boltzmann destribution. 

The variation of population density as a function of level 
energy is quite different from what is usually observed in the 
cold ISM. Typically the first two rotational levels are consis- 
tent with a rather steep slope with a temperature 7b i ~ 80 
K while the higher rotational levels flatten, giving the ap- 
pearance of a higher "temperature" of several hundred de- 
grees or more (ISpitzer & Cochra3 l973j). The excess column 
in the J > 1 levels has long been thought to be caused by 
far-ult raviolet continuum fluoresence, but recently iGry et all 
(2002) has questioned this hypothesis. They favor a true two 
temperature model for the ISM. 

Regardless, here we have the opposite case. In each vi- 
brational level v" = 0, 1 the trend is for the lower rota- 
tional states to have flatter slopes (higher temperature) than 
the higher rotational states. Reproduction of the energy 
level distribution found here will be an interesting chal- 
lenge for H?_ excita tion models (c.f. Smtzer & Zweibel ll 19741 
Black & Daigarnolll976t IShulllll978t Ivan Dishoeck & Blackl 
19861: ISternberg & Dalgarnoll 19891: iDraine & Bertoldilll996l) . 



We note how the modeling by Spit zer & Zweibell (1 19741) 
shows qualitatively how the "curvature" of the population 
density can go from "concave" to "convex" as the parame- 
ters change from a low temperture, density and photoexcita- 
tion rate environment to either a high photoexcitation rate or 
a high density environment. However, they do not attempt to 
explore temperatures in excess of 1000 K. 

4.2. Central Star Spectral Energy Distribution and 
Extinction 

SED models with CS parameters taken from the literature 
provide a fairly precise match to the stellar continuum ob- 
served by FUSE. This allows us to determine the line-of- 
sight extinction, and provides a reasonable stellar continuum 
to use for assessing the success of our model of atomic and 
H2 absorption. The combined SED and ab sorption model 
has been used by iMcCandliss & Krukl d2007l) to identify 63 
photospheric, nebular and non-nebular absorption features of 
ionized and neutral metals, lurking amid the sea of hydrogen 
features. 

In Figure @] we show a log-log plot of the CS SED from 
900 - 20000 A, as measured using FUS E IUE, and the APO 
PIS, along with the optical photometry dTylenda et al1ll99 lb 
ICiardullo etTniigggtlBenedictetail^OOSh . Over-pl otted is 
a synt hetic stellar flux interpolated from the grid of iRauchl 
(2003) with logg = 6.5, T = 120,000 K, and ratio of H/He = 
10/3 by mass. This model includes no metals and is consis- 
tent with, a lthough sligh t ly hot ter than, the quantitative spec- 
troscopy of Napiwotzki ( 1999). It is slightly cooler than the 
determination of iTraulsen et alJ (120051) . The temperature and 
gravity adopted for our model has slightly less pressure broad- 
ening and is a better match to the observed Lyman lines to- 
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wards the series limit. We also adopt a stellar mass of 0.56 
M© as suggested by post asymptotic giant branch evolution- 
ary tracks. 

Use of this mass along with the above gravity required a 
distance of 466 pc to match the absolute fl ux, which is at the 
upper limit given by Be nedict et al.1 (120031) . Beyond our im- 
mediate need to match the absolute flux for the given gravity 
and mass there is n o particular reas on to prefer this distance 
over that derived by Ben edict et al.1 Questions regarding the 
acceptable uncertainty in distance, absolute flux and derived 
stellar parameters are best left for a stellar model specifically 
tailored to include the effects of metals, gravity and evolu- 
tionary state. For our purpose, the high stellar temperature 
places the SED longward of the Lyman limit in the Rayleigh- 
Jeans regime, so the shape of the SED is insensitive to our as- 
sumptions of temperature, gravity, metallicity, mass and dis- 
tance. Consequently, long-range changes in SED shape pro- 
duced by reddening can be constrained with a high degree of 
confidence. 

The model flux is shown in Figure |4] with reddenings of 
E( B-V) = 0.00 0.005 and 0.05 assuming the Ry = 3.1 curve 
of ICardelli etall d 19891) . An extinction of E(B-V) = 0.05 is 
clearly too high. As a convenience we adopt a very modest ex- 
tinction of E(B — V) = 0.005, as seen most clearly in Figure Q] 
to make the continuum towards the Lyman edge match the 
stellar SED. The lack of any significant extinction is rou ghly 
consistent with the total nebular hydrogen, derived in § 14.31 
and the standard conversion for color excess to total hydro- 
gen column of N(H )/E(B-V) = 5.8 x 10 21 cm~ 2 mag _1 
ilBohlin et al.lfl978h . 

4.2. 1 . Balmer Line Reddening 

Reddening by dust in PNe is often estimated by compar- 
ing the measured Balmer emission ratios to the intrinsic ratios 
produced by hydrogen recombination. The intrinsic intensity 
ratio (in flux units) of 7(//a)//(Hp) ~ 2.859 at an electron 
temperature of 10 4 K and density of 10 2 cm~ 3 , is accurate to 
within « ± 5% for a wide range of the electron temperature 
and density, (5000 < T e (K) < 20000, 10 2 < « f (cm- 3 ) < 10 6 ) 
dBrocklehurstll 1 97 Tl) . The extinction parameter c is given in 

I(k)/I(H$) =F(X)/F(H$)l0 cm (1) 

where the intrinsic ratio is I(X)/I(H$), the observed ratio 
is F(X)/F(H$) and f(k) is the line attenuation, relative to 
HP, at any wavelength for the give n extinction curve (c.f. 
iMiller & Mathewslll972T: ICahnll 19761) . The extinction param- 
eter may be converted to selective exti nction by adopt i ng the 
ratio of c/E(B-V) = 1.5 as shown bv lCiardullo et al l (1 19991) 
in their Figure 5. 

Balmer emission line spatial profiles for Ha, Hp, Hy, and 
H8 in a ± 140" region to either side the CS at a position angle 
of 35 °, were acquired as discussed in § 13.31 The line ratios, 
with respect to are displayed in Figure [5] with a two pixel 
(0"9) binning. The reddenings derived from these ratios (Fig- 
ure [5]l are inconsistent and are summarized in Table [3] The 
first column specifies the line ratio, the second is the extinc- 
tion multipler taken from iBarkerl (|l984), the third is the in- 
trinsic line ratio for a te mperature of 10,000 K and a density 
of 100 cm~ 3 taken from Bro cklehurstl (119711) . the fourth and 
fifth columns are the average line ratios and standard devia- 
tions for the SW and NE regions after rebinning the data by 
sixteen pixels (7 '.'2), and the sixth and seventh columns are 
the derived selective extinctions and associated errors. The 
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See figure[7j ' Emission columns are upper limits. Error is ± 0. 15 



dex.* Absorption columns used in H I model. 

NE Ha/Hp ratio suggests E(B-V) = 0.10 ± 0.02, while all 
the rest of the ratios are consistent with much lower or no ex- 
tinction. 

The line ratio averages in the fiduical SW and NE regions 
have standard deviations of 1.5 - 3 % after rebinning by six- 
teen pixels. The large error bars in the derived E(B— V) serve 
to emphasize that Balmer line ratio determinations to a pre- 
cision of better than 1% are required for E(B— V) < 0.01. 
We note that in the immediate vicinity of the star Ha/HP = 
3.03 yielding E(B— V) = 0.051, which is much too large to be 
consistent with the observed far ultraviolet SED (see § 14.21 ). 
assuming the line-of-sight extinction is close to the Galactic 
standard. 

We concluded in the previous section that dust along the 
diffuse line of sight produces E(B — V) < 0.01 mag. Here we 
find the extinction in the extended medium is low to nonex- 
istent. We conclude that dust is not widespread, although it 
may be confined to the clumpy medium. The inconsistancies 
in Balmer reddening measures suggest that a process other 
than extinction by dust is contributing to Ha/Hp ratio excess. 
We will discuss this issue further in § 15.41 

4.3. Atomic Hydrogen Absorption Model 

In Table |2] we list the column densities, doppler velocities 
and heliocentric velocity offsets inferred from the Dwingeloo 
data (§ 13.4b . The 21-cm emission column densities listed 
in column 2 have contributions from both background and 
foreground emission, so they provide an upper limit to the 
expected foreground absorption column densities. We are 
mostly concerned with finding reasonable numbers for the 
nebular absorption, i.e. those components shortward of the 
stellar systemic velocity of —42 km s ; components 1-4. 
We will refer to components 5 - 9 as the non-nebular ISM 
components. 

The adopted neutral hydrogen absorption model, the com- 
ponents of which are listed in column 5, was determined with 
the following procedure. We initially allowed the two compo- 
nents nearest to V sys (5 and 6) to retain the maximum column 
suggested by the emission line fit. The two components nicely 
filled the centers of the Lyoc - Ly/j lines. We then moved to the 
blueward side and reduced the column of component 1 until a 
good agreement was found in the most saturated lines Lyr| - 
Ly/j. We continued to add in components 2, 3 and 4 reducing 
the columns as necessary to smoothly match the absorption in 
the saturated lines. Once we were satisfied with the blue edge 
fit we proceeded to the redward components. It was deter- 
mined that component 9 produced too much red edge absorp- 
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FIG. 10a. — Lyman series model fits (Lya - fi). The long dashed lines near the center of the profiles mark the velocities of the absorption components 
determined from the Dwingeloo data. In Lya the two STIS E140M orders are shown in orange and green. The FUSE data are shown as orange (sl2) and 
green (s!2). The stellar c ontinuum is in purple and the absorption model is in black. The numbered +'s mark absorption features tabulated in Table 2 of 
McCandliss & Kruk (2003). 
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FIG. 10b. — Lyman series model fits (LyT| -/j). Colors and marks are the same as in Figure [TOa! 
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TABLE 3 

Extinctions Derived from Balmer Ratios * 





7W 


Int Ratio : 


Obs Ratio SW* 


Obs Ratio NE :|: 


E(B-V) SW 


E(B-V) NE 


Ha 

w 

Hy 
HP 
HS 
HP 


-0.33 
0.15 
0.20 


2.859 
0.4685 
0.2591 


2.868 ± 0.041 
0.4627 ± 0.0093 
0.2585 ± 0.0075 


3.205 ± 0.066 
0.4724± 0.0144 
0.2508± 0.0082 


0.003 ± 0.013 
0.024 ± 0.039 
0.003 ± 0.042 


0.100 ±0.018 
-0.016± 0.059 
0.047 ± 0.048 



* Data rebinned by 16 pixels* f(X) from Barker 11984).' Instrinsic line ratios from Brocklehurst 
fi97lh .* SW values are averaged over (—80,-50)", the NE over (50,80)"; see dashed lines in Figure[5] 



tion in all the Lyman series lines and it was eliminated from 
further consideration. We found that only a modest amount 
of component 8 was necessary to match the redward edge of 
Lya. Finally component 7 was added to reduce the absorp- 
tion gap between components 6 at -20 km s _1 and 8 at 11 
km s . Changing the column density of any one component 
by the error in the emission line fit (0.15 dex) did not appre- 
ciably cause deviations of the resulting model with the data. 
The constraints on the formation and destruction of molecu- 
lar hydrogen within the nebula, which we discuss in § 15.31 
are immune to the details of the velocity distribution of H I 
derived f rom t he Dwingeloo data. 

Figure UOal shows an overlay of the adopted atomic and 
H2 absorption model on the Lya - Ly/j lines from STIS and 
FUSE. The Lya profile from the STIS E140M observation 
spans two orders, plotted as orange and green. The rest of the 
Lyman series are from FUSE with the sl2 spectrum shown in 
orange and s21 shown in green. The stellar SED is shown in 
purple and the model is in black. The locations of the the indi- 
vidual H I velocity components from Table [2] are marked with 
vertical dotted l ines. The model includes the H2 absorption 
determined in § 14.11 A comparison of the hydrogen absorp- 
tion model against the sl2 and s21 spectra for the entire FUSE 
wavelength range along with metallic absorption system iden- 
tifications, arising from the photosphere, the nebula and non- 
nebular ism, can be found in lMcCandliss & Krukl(l2007l) . We 
see that the model continuum in the vicinity of Ly(3 is ~ 10% 
too high, while at the shortest wavelengths it is low by a sim- 
ilar amount. Ther e are also indications of ab sorption from O I 
and N II features (iMcCandliss & Krukll2007l Table 2). 

Although the agreement of the model with the data is excel- 
lent, the procedure adopted here for constraining the neutral 
hydrogen column in the nebula is not ideal and in all likeli- 
hood provides a non-unique solution. However, the result is 
at least plausible, well within factors of two, given the uncer- 
tainty associated with the large background subtraction ne- 
cessitated by the coarse angular resolution of the Dwingeloo 
data cube. A high angular resolution 21 cm mapping of the 
nebula would provide a more useful constraint on the velocity 
components, column densities and doppler widths associated 
with the nebula. We are mainly concerned with the total H I 
column density of the nebula rather than the individual com- 
ponents. The sum of components 1 - 4 is N(H \) ne b = 3.0 ± 
1 x 10 19 cm 2 . The non-nebular components 5-8 total to 
N(H l) non = 6.2 ± 1 x 10 19 cm 2 . 

4.4. CO Upper Limit and the C I Pressure Diagnostic 
AAA. CO Upper Limits 

iBachiller et al.1 (12000)) presented a detailed contour map of 
the CO (2-1) emission in M27. The contours fall steeply to- 
wards the CS where, using the conversion factor suggested 



by iHuggins et al.1 0996), we estimate an upper limit for the 
CO column of 10 15 cm~ 2 along the line-of-sight. CO has 
a number of strong absorption bands in the far UV, such as 
the numerous A-X bands spread throughout the STIS E140M 
bandpass, and this column density would easily be detected. 
In Figure QT| we show the E140M spectral region surround- 
ing the A-X (2-0) 1477.565 A bandhead. In the top panel 
we overplot CO absorption models for representative excita- 
tion temperatures of 3000 K and 300 K at a column density 
of 10 15 cm~ 2 in blue and red respectively. We assumed the 
doppler parameter is the same as for H2, 6.5 km s . Clearly 
the column density is much lower than 10 15 cm~ 2 . In the 
lower panel we show models for 3000 K with a column den- 
sity of 4 x 10 14 cm~ 2 and for 300 K with a column density of 
8 x 10 13 cm~ 2 . The structure in the models is on the order of 
the noise in the observation, so we will use these values as the 
(excitation temperature dependent) upper limits for line-of- 
sight CO. Using a lower excitation temperature would yield 
an even lower the upper limit. 

4.4.2. C / Fine Structure Pressure Diagnostic 

Although CO was not detected we do detect excited C I, 
C I* and CI** multiplets in the E140M data with a velocity 
of -75 km s , coincident with the velocity of the excited H2. 
The level populations of the three fine structure levels in the 
ground of C I, assuming a long enough time has passed for 
equilibrium to be established, is a detailed balance of the de- 
excitation and excitation rates set by collisions of C I with 
other particles (the impactors - H I, He I, p, e~, ortho-H2, 
para-H 2) and the radiative dec ay rates of the fine structure 
levels. iJenkins & Shayal (1 19791) present a useful diagnositic 
for gas pressure based on the ratios of the column densities 
for Nci* and Nci** relative to the total column density N to t= 
Nci + N a * + N a **, (fl = N C iJN M and f2 = N C i**/N tot ). The 
pressure is constrained by comparing the measured values of 
fl and f2 against theoretical loci for fl and /2 at variable 
temperature and density. 

We determined fl and f2 by % 2 fits to the absorption pro- 
files of the C I A11656.27 - 1658.12 UV2 multiplet using the 
iMortonl d2003l) oscillator strengths. The result is shown in 
Figure [12] We determined a doppler parameter of b = 4.5 
km s- 1 and we find N a = 3.0 ±0.5 x 10 13 cirr 2 , N a * = 
8.5 ±0.5 x 10 13 cm" 2 and JV C /« = 10.5 ±0.5 x 10 13 cirr 2 , 
yielding (/l,/2) = (0.386, 0.477) ± 0.04. 

The theoretical l oci for fl and f2 depe nd on the mix of 
impactors assumed. IJenkins & Trippl (1200 ll) discuss three ex- 
tremes: Casel where all the impactors are neutral hydrogen, 
Case 2 where all the hydrogen is molecular and Case 3 where 
all the hydrogen is ionized. They find Case 2 and Case 1 
to yield very similar loci in the (/l,/2) plane. Using their 
code (kindly provided by Jenkins, private communication) we 
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FIG. 1 1 . — HSTSTIS E140M spectrum of the region around CO A-X (2-0), the strongest band in the A-X absorption series. Top - Two different excitation 
temperature models (3000 K, blue) and (300 K,red) of the (2-0) band using the same column, 10 15 cm -2 , as suggested from CO mm measurements. Bottom - 
noise bounded upper limits to the line-of-sight CO column for the two excitation temperatures. Note the scale change between the two panels. The upper limit 
drops for decreasing excitation temperatures, because the population becomes concentrated in the lower rotational levels as temperature decreases. 



show the results for Case 1 (left) and Case 3 (right) in Fig- 
ure [13] The crosshairs mark the (/l,/2) point and the figure 
ranges over the ± 0.04 error limits. Crosses mark the density 
at intervals of 0.1 dex. The log of the pressure, in units of 
log(n/ t T) = log(P/k), is printed for every density point along 
the loci of points for each constant temperature. 

For the neutral case the (0.386, 0.477) point lies on the 
log(T) = 2.6 curve with a pressure of log«/,r = 5.9, while 
the ionized case has log(T) w 3.0 with a pressure of 4.9 < 
logn/,r < 5.0. The allowed range of pressure for the ± 0.04 
error box in the neutral case is 5.7 < log n/,r < °° and for the 
ionized case 4.8 < logn&r < °° with the °° density case being 
associated with the coldest allowable logarithmic temperature 
of«2.1dex. 7 

The coincidence of the C I velocity with the H2 velocity 
argues that these two species are co-located. The ~ 2000 K 
implied by the H2 ro-vibrational distribution is closer to the 
temperature given by the ionized case than the neutral case. 
This suggests the excited C I and H2 are located in a warm 
and electron rich medium with n e ~ 80 cm~ 3 . In this envi- 
ronment it may be possible for H2 to form at the interface 
between the ionized medium and the PDR clumps by H + e~ 
— > H + hv followed by H + H~~ — > H2 + e~ mechanism , 
(c.f. lAleman & Gruenwaldll2004l: iNatta & Hollenbach|[l998l) . 
We note that our total column density for the neutral carbon, 
Ncitot = 2.2 x 10 14 cm~ 2 is ~ 2 orders of magnitude higher 
than what is expected for Ncitot in typical H II reg ions, as dis- 
cussed in Appendix A of Jen kins & Shaval (119791) . It is likely 
that the source of C I in the diffuse medium is CO dissociation 

7 For constant temperture curve the (/l,/2) points converge to a point as 
the logarithmic density approaches 6 and the diagnostic becomes less dis- 
criminatory. 



in the clumpy medium. 

In addition to UV2, we also examined the fits for the UV3, 
UV4 and UV5 multiplets. We found that the % 2 fits to these 
multiplets yielded slightly different results. In the cases of 
UV3 and UV5, there were intervening unidentified absorp- 
tion features that caused the % 2 to drift to higher columns for 
the Nci** lines. In the case of UV4 the column densities were 
different but the fl and /2 ratios were essentially t he same. 
Currently there are disagreem ents in the literature (c.f. lMortonl 
l2003t iJenkins & Trippll2001h regarding the relative strengths 
of the higher multiplets with respect to UV2. Since the result- 
ing flandfl ratios all lie within the ± 0.04 error box, we de- 
cided to rely upon the multiplet with the most certain oscilla- 
tor strengths, UV2. Our conclusions that the C I is excited, the 
column density is high compared to typical H II regions, and 
the temperature implied by the ionized impactors are closer to 
that implied by the H2 level population, are unchanged by our 
reliance on UV2 alone. 

4.5. FUSE and STIS Absorption Profiles - Ionization 
Stratification 

The metal line profiles show absorption from velocity com- 
ponents that we associate with the stellar photosphere located 
near the systemic velocity of the star, the blueshifted nebu- 
lar outflow, and the redshifted non-nebular ISM. In general, 
the high ionization species tend to appear with photospheric 
and nebular velocity components while the low ionization and 
neutral species tend to show nebular and non-nebular compo- 
nents. A gravitational redshift of « 5 km s _1 can be seen in 
the high excitation stellar O Ve A-1371 line. (We use the "e" 
specification with the atomic species designator to indicate 
transitions to lower levels other than the ground state). 

We show a select subset of the metal line absorption profiles 
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in Figures[l4]-[T7] FUSE sl2 and s21 spectra are shown in 
orange and green respectively, while the STIS lines are plotted 
in purple (and grey if two orders contain the same wavelength 
range). The continuum model, including the hydrogen lines, 
is plotted as a thin black line. The heliocentric reference frame 
is used. The black dashed vertical lines mark the location of 
the hot H2 component at -75 km s the systemic velocity at 
-42 km s _1 , and the cold H2 component at -28 km s . The 
red dashed vertical line marks the gravitational redshift of the 



systemic velocity. To convert from the nebular rest frame, 
subtract the heliocentric velocity from the systemic velocity. 

4.5. 1. Stellar Photospheric Features and the Absolute Wavelength 

Scale 

Comparison of the FUSE and STIS wavelength scales re- 
vealed a systematic offset. The reconciliation of this offset is 
essential for investigating the kinematics of the nebular out- 
flow, where we seek to determine the velocity of the various 
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molecular and atomic features with respect to the systemic 
velocity of the nebula. Lines that arise from the photosphere 
should match the systemic velocity of the system less the off- 
set caused by gravitational redshift. An absolute reference 
to the heliocentric velocity was established by close exami- 
nation of the O Ve XI 37 1.296 feature, shown in the middle 
panel of Figure [14] This narrow line results from a transition 
between two highly excited states in O Ve and is expected to 
be an excellent indicator of the photospheric restframe (Pierre 
Chayer private communication). For a compact object of this 
mass and radius (§|2]i we expect the photospheric lines to ex- 
perience a gravitational redshift, V gr = c((l — 2 ™i ) ~ % — 1) = 

5.1 km s . Applying a shift of-13 km s _1 to the STIS spec- 
trum placed the centroid of the O V A. 1 37 1.296 at - 37 km s~ 1 , 
as expected for a V sys of —42 km s -1 ([Wilson 1953]). 

Examination of overlapping FUSE and STIS spectra in the 
wavelength regions below 1190 A revealed some systematic 
disagreem ents. In the origin al analysis of the FUSE M27 
spectra by McCandliss (2001) the hot nebular H2 component 
was defined to be at -69 km s _1 . We found it necessary to 
shift the FUSE spectra blueward by -6 km s , such that the 
hot nebular H2 is now at -75 km s _1 . The most useful overlap 
lines for assessing the alignment were the doublet blend of 
C IVe Ul 168.849, 1168.993, bottom of Figure [H and the 
the narrow O Vie Ul 171.56, 1172.44, top of Figure [H We 
note the wavelengths of O Vie doublet given in the National 
Institute of Standards and Technology (NIST) online tables 
( |http: //physics .nist .gov/Phys Ref Pat a/ ASP/ lines-form.html) 
appear to be in error by sa -0.42 A fsee Uahn et alJl2006h . 
These spectra are comparatively noisy, but the alignment with 
the FUSE spectra agrees as well as the alignment between the 
sl2 and s21 spectra. We also used the FUSE N I multiplets at 
Ul 134 - 1 135 and the STIS N I multiplets at U1200 - 1201 
along with excited H2 lines that appear in the STIS bandpass 
ab ove 1190 A with the con tinuum plus hydrogen model 
of iMcC andliss & Krukl ((2007) to assess the consistency of 
the wavelength reconciliation. The error in the systemic 
velocity (§ EJ is of order the FUSE resoution element and 
is three times the STIS resolution element. We consider the 
agreement of line profiles from spectra acquired with two 
different instruments to be excellent. 

4.5.2. Photospheric + Nebular Features 

O VI, N V and C IV — Figure Q3] shows the high ioniza- 
tion resonance doublets of O VI U1032.62, 1037.62, N V 
U1238.821, 1242.804 and C IV U1548.204, 1550.781. All 
these lines show photospheric absorption to the red of the sys- 
temic velocity. The strength and terminal velocity of the blue 
shifted portion of the line profiles increases with decreasing 
ionization potential. The O VI lines show only slight signs of 
blue shifted nebular absorption. The nebular absorption com- 
ponent in the N V lines is strong only between -42 and -75 
km s just reaching saturation at ss -60 km s . In contrast, 
the nebular absorption component in the C IV lines spans —42 
to -1 15 km s _1 and is completely saturated from -50 to -95 
km s . 

4.5.3. Intermediate Ionic and Neutral Nebular Absorption 

C III - II, NIII- I and 01— Figure [16] show lines of C m - 
H, N III - I and O I. Like the C IV lines, the C III - II lines 
are heavily saturated throughout the nebular flow region blue- 
ward of -42 km s _1 . The saturation makes it difficult to tell 
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FIG. 17. — Comparision of the Si IV - 2 to S IV - 2. The low ionization species, 
high i oniz ation species, Si IV and S IV, dominate near -60 km s _1 , Transition in 
Figure [Ml fbr a description of the colors. 



Si II and S II, dominate in the zone between -75 and -1 10 km s 1 while the 
ionization occurs where absorption by the neutrals and H2 is strongest. See 
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whether any one ion is dominant in the nebular outflow. N III 
is strongly blended with overlapping H2 features. However, 
it shows saturated absorption between -60 and -100 km s . 
N II A.1083.994 shows absorption throughout the flow, being 
less saturated at low velocities and becoming completely sat- 
urated at -75 km s _1 . This line also shows stronger nebu- 
lar absorption than non-nebular absorption. In contrast, N I 
A,1200.223 shows weaker nebular absorption than non-nebular 
absorption. Both N I A,1200.223 and O I show an absorption 
component centered on -75 km s~ 1 , which coincides with the 
velocity of the absorption maximum in H2 and C I. 

57 IV - II, and S IV - II — The transition zone between high 
ionization and low ionization occurs at the velocity of -75 
km s _1 where H I, C I, N I, O I and H2 show up most strongly 
in the nebula. The high ionization - low velocity, low ioniza- 
tion - high velocity dichotomy is best illustrated by examin- 
ing the velocity profiles of the low abundance metals Si and S 
in Figure [17] Si II is stronger than Si IV in the zone between 
-75 and -1 10 km s _1 , while Si IV is stronger than Si II in the 
zone between -42 and -75 kms" 1 . The Si III A.1206.500 line 
is saturated throughout most of the flow. The S IV - S II ions 
show similar behavior. 

5. DISCUSSION 
Here we enumerate the main findings from the analysis. 

1. The nebular H2 is highly excited with a ground state 
ro-vibration population that deviates significantly from 
the best fit single temperature Boltzman distribution of 
2040 K. The total column density is 7.9 x 10 16 cm" 2 . 
The deviations are characterized by a flatter slope for 
the lower rotational levels (0 < J" < 5) than for the 
higher rotational levels (6 < J" < 11). 

2. Continuum fluorescence of H2 has not been detected to 
the background limit of the sounding rocket observa- 
tion w 5 x 10~ 17 ergs cm~ 2 s _1 A -1 arcse c~ 2 . How- 
ever, L yq fluorescence has been detected by Lup u et all 
(2006). This emission mechanism is possible only for 
H2 with a significant population in the v" = 2 vibrational 
state in the presence of a strong Lyoc radiation field. 

3. The stellar SED exhibits negligible stellar reddening 
along the line-of-sight. This finding is at odds with the 
variable reddening indicated by the ratio of Ha/Hp in- 
tensity as reported in the literature and confirmed in this 
study in the immediate vicinity of the star. We note that 
the reddening determined by Hy/H(3 and H8/H(3 is con- 
sistent with zero and conclude it is possible that another 
physical process, aside from reddening by dust, is caus- 
ing the discrepant Hoc/Hp. 

4. The nebular H2 appears at a heliocentric velocity of - 
75 km s , corresponding to an expansion velocity of 
33 km s _1 in the nebular restframe. H I, C I, N I, and 
O I exhibit strong resonance line absorption at the same 
expansion velocity. The molecules and neutrals de- 
marcate a transition velocity between a regime of high 
ionization and low expansion velocity (10 < v/„ < 33 
km s ) and a regime of low ionization and high ve- 
locity (33 < \'i ow < 65 km s ). The upper limit to 
the terminal expansion velocity is as 70 km s _1 . These 
observations, along with the absence of any signs for a 
high velocity radiation driven wind, are a challenge to 



interpret in the context o f the simple interac t ing winds 
model for PN fo rmation dKwok et al.lll978l). IMe aburn 
d2005l) (see also IMeaburn et aT[ T2005a; Meaburn et al. 
2005b), who used high resolution position-and-velocity 
spectroscopy of the optical emission lines to derive the 
ionization kinematics of several objects, have arrived 
at a similar conclusion. They find that ballistic ejec- 
tion could have been more important than interacting 
winds in shaping the dynamics of PNe. Bali ck & Frankl 
(2002) have reviewed in detail some of the more so- 
phisticated thinking about the dynamical processes that 
shape PNe. 

5. Analysis of the C I multiplet fine structure populations 
gives better agreement with the expected temperature of 
~ 2000 K for collisions with protons and electrons as 
opposed to collisions with neutrals or molecules. This 
suggests the molecular and neutral transition compo- 
nent is in close contact with ionized material. 

6. The upper limits to CO in the diffuse line-of-sight 
gas range over 8 x 10 13 cm" 2 < N(CO) < 4 x 10 14 
cm~ 2 depending on the assumed CO rotational excita- 
tion temperature. The high column density of C I found 
in the nebula is consistent with it being a daughter prod- 
uct of CO dissociation. 

7. The lower limit to the neutral to H2 column density ratio 
is 7V(H Y)/N(H2) = 127 in the transition region (using 
component number 3 in Table |2). 

We will now discuss the constraints imposed by these find- 
ings on the ionization kinematics, mass structures, nebular 
dust, and the excitation, formation and destruction of H2 
within the nebula. 

5.1. Ionization Kinematics 

Low mass to intermediate mass stars begin to evolve into 
AGB stars when core nucleosynthesis ceases and a degenerate 
core forms. The AGB phase lasts « 250,000 years and is 
marked by a series of high mass loss events that merge to 
produce a slow wind. At some point the star collapses into a 
hot compact degenerate object. In the PN formation scenario 
suggested by lKwok etall (119781) . a fast (~ 1000 km s _1 ) low 
density radiation driven wind from the hot star shocks and 
ionizes the slow moving (~ 10 km s ), high density mostly 
molecular AGB wind, resulting in the expanding structures 
w e see at latter times. 

iNatta & Hollenbachl 0998) have computed detailed time 
dependent changes in the atomic and molecular emissions in a 
phenomenologically motivated interacting wind model. Their 
model consists of an interior ionization region bounded by a 
shell propagating outward with a velocity of 25 km s over- 
running and shocking a mostly molecular AGB wind with a 
velocity of 8 km s _1 . An ionization-bounded time-dependent 
PDR results with three layers. These are, a fast moving ion- 
ized shock interface, an adjacent dissociation layer, and a slow 
moving molecular layer further downstream. They find most 
of the shell mass to be vibrationally excited H2 and H I pro- 
duction to be modest. After w 4000 years they find the in- 
frared H2 line diagnostic of S(l) (1-0) / S(l) (2-1) indicates a 
switch from thermal excitation to fluorescent excitation. They 
also find that H2 does not co-exist with CO, as dissociation fa- 
vors C + + O. Consequently, they suggest molecular masses 
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calculated for old PN based on CO emission and standard 
CO/H2 ratios may be underestimated. 

M2 7 has a kinematic age of « 10,000 yrs dO'Dell et alj 
2003), so perhaps it is reasonable that the nebula we see today 
bears little resemblance to this simple layered, two-velocity 
interacting wind picture. We see no evidence for a shell of 
predominately molecular material along the diffuse line-of- 
sight medium, nor do we see neutral or molecular material at 
low expansion velocity as would be expected for an unper- 
turbed AGB wind. The upper limit to the molecular fraction 8 
in the neutral velocity component expanding at 33 km s~'is 
/h 2 = 0.016, assuming N(H)=N(H I). We do see evidence for 
velocity structure in the intermediate ionization states (Fig- 
ure s[T6"l-[T7l) and H I (Tabled out to expansion velocities of 
w 65 km s -1 , but there is no evidence for P-Cygni profiles 
with ~ 1000 km s terminal velocities in any of the high 
ionization lines (Figure \T5[. 

We find instead a rather permuted velocity structure of neu- 
tral and molecular material embedded (in velocity space) be- 
tween the high and low ionization medium. Our findings sug- 
gest that the AGB wind has been accelerated through interac- 
tion with the hot star and that much of the original mostly 
molecular atmosphere has been almost completely dissoci- 
ated. The surviving molecular material is currently interact- 
ing strongly with the highly ionized medium and we postulate 
that the intermediate ionization species are created and accel- 
e rated in this interaction . 

IVillaver et al.l (l2002bl) have treated the radiative hydrody- 
namical problem with great rigor for a range of progenitor 
masses, using as a starting point for their PN si mulation the 
structur e found at end of AGB evolution by IVillaver et alJ 
(2002a). They find a main nebular shell develops surrounded 
by a relatively unperturbed halo. They give detailed figures 
for the evolution of density, velocity and Ha emission bright- 
ness. We note that the amplitude of the velocity structures 
for their 1-2.5 M sun models are in the range of the out- 
flow velocities exhibited by the line profiles in Figures [TBI - 
[TTl Unfortunately they give no details on the ionization kine- 
m atics of the main and halo regions as a function of time, 
as iNatta & Hollenbachl (Q998) have done. A model com- 
bining the detailed time de pendent calculation of the atomic 
and molecular emissions of INatta & Hollenbachl dl998|) with 
the radiative hydrodynamical rigor of IVillaver et al.l ([2002b) 
would be useful for interpreting these observations. 

5.2. Mass Structures 

Identifying kinematic components appearing in spectra 
with spatial structures in images must be approached with 
care. The computed total mass of these velocity components 
will depend on the assumed diameter and thickness of partic- 
ular structures in the nebula. However, it seems reasonable 
to expect the 33 km s _1 transition zone, where most of the 
molecular and much of the neutral hydrogen appears, to be 
at or near the region where the Ha brightness drops to near 
zero. Inside this zone we expect the medium to be com- 
pletely ionized. The minor axis of the bright Ha in M27 is 
« 300" and the full diam e ter in cluding the faint halo detected 
by Papama storakis et ail (1 19931) is 1025". Assuming a dis- 
tance of 466 pc gives radii for the transition region and the 
total nebular extent to be w 0.34 and 1.16 pc respectively. 
These radii agree reasonably well with the bright Ha emit- 

8 By nucleon number /g 2 = 2N(H 2 )/(2N(H 2 ) + N(H)) with 
N(H)=N(H I)+jV(H II). 



ting region and the full extent of the AGB wind found in the 
IVillaver et all (l2002albl) simulations after w 10,000 years . 

We estimate the masses of the various structures given the 
following (gross) assumptions. Let the ionized region be a 
sphere with a radius of 0.34 pc and an electron density of ~ 
300 ± 200 cm~ 3 (Barker 1984). The ionized mass becomes 
0.4 < Mion < 1-9 M©. If the total neutral hydrogen column 
density (3 x 10 19 cm~ 2 ) is uniformly distributed in a shell 
between 0.34 and 1.16 pc with a mean density of 12 cm~ 3 , 
then M s h = 1.7 M . Our constraint on the progenitor mass, 
the sum of the ionized, shell and degenerate star masses, is 
2.6 <M pro <4.2M Q . 

By comparison, the molecular material in the diffuse neb- 
ular medium has a much lower mass. Although the H2 is 
associated with only one neutral hydrogen velocity compo- 
nent, we assume it is spread uniformly throughout the 0.82 
pc thick shell. Using this maximal volume yields an upper 
limit on the molecular mass in the diffuse medium, Mdmol S 
(1.7 M Q )/375 = 4 x 10~ 3 M@. iMeaburn & Lopezl ( fl99l 
have estimated the total mass of the eleven largest globules, 
which appear as dark knots against the bright nebular O III 
emission, to be M omo \ = 1.1 x 10~ 3 M Q . They used mea- 
surements of O III transmission and assumed the general ISM 
gas-to-dust ratio dBohlin et al.l 1 1978h . The Meaburn et al. 
clump mass e s timate s compare reasonably well to those of 
Huggins et al. ( 1996) who estimate the total molecular mass 
associated with the CO emitting clumps to be M am ,i = 2.5 x 
10~ 3 M Q . They assumed a thin ring geometry and a fraction 
of CO/H 2 = 3 x 10" 4 , a valu e > 10 times than the high- 
est found by lBurgh et all (120071) for the diffuse ISM. Use of a 
lo wer ratio will increase the CO clump mass estimate as noted 
bv lNatta & Hollenbachl <fl998h . 

Our upper limit to the diffuse molecular mass is on the same 
order as the clump mass estimates. The upper limit on the to- 
tal molecular mass in the nebula, the sum of the diffuse and 
CO derived masses (assuming the high CO/H2 ratio), is M,„„/ 
= 6.5 x 10~ 3 M . We conclude that there is no large mass of 
molecular material in M27, eithe r in the diffuse or clu mped 
medium. This is in contrast to the lMeixner et alj d2005l) find- 
ing for the Helix, for which they estimate the mass of all the 
H2 emitting globules to be ~ 0.35 M@. We note that all of 
the molecular mass estimates for M27 are on the order of the 
mass of the planets in our own solar system, w 1.3 X 10~ 3 
Mq. The suggestion that dense circumstellar bodies have sur- 
vived the AGB phase and may constitute a significant fraction 
of the molecular material in the nebula environment cannot be 
ruled out. 

The absence o f CO absorption in th e diffus e medium is con- 
sistent with the INatta & Hollenbachl d!998l) finding that CO 
and H2 are unlikely to co-exist in the gas phase, because CO 
is quickly diss o ciated after leaving the vicinity of the clumps. 
Huggins et al. (2002) confirmed this picture with exquisite 
observations of the globules of the Helix. They show how 
CO survives only in the shadows of the clumps. The heads of 
these globules show strong coincident emission of H2 S(l) (1- 
0) and Ha, and we a ker co incident emission in the shadowed 
zone. lO'Delletail d2003l) note the globules in M27 tend to- 
wards a more irregular and tail-less appearence than those in 
the Helix. Nevertheless the observation of H2 emission in the 
vicinity of the globules suggests they may be a source of H2 
for the diffuse nebular medium. 

The temperature of the diffuse medium inferred from the 
molecular hydgrogen ro-vibration levels exceeds the subli- 
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matio n temperature of refractory dust, ~ 1500 K (Glassgold 
1996). While it is unlikely that the gas and dust will be in 
thermodynamic equilibrium, the absence of any significant 
reddening of the central star SED suggests that the diffuse 
environment is as inhospitable to dust as it is to CO and H2. 
Dust, or perhaps even compact bodies, may e xist in the CO 
rich dense clumps, which we will argue in § 15.3.21 are the 
likely source of H2 for the diffuse medium. 



5.3. 



Constraints on the Excitation, Formation and 
Destruction of H2 

5.3.1. H2 Excitation Processes 



The models of iNatta & Hollenbachl (1 19981) suggest that in 
the latter stages of planetary nebula evolution the H2 is being 
excited via far-UV continuum fluorescence, yet we have failed 
to detect this emission. We briefly review this process. 

H2 is excited from the ground vibrational level of X 

into higher electronic states, predominately B and C x Yl u 
(the Lyman and Werner bands), by the absorption of a far- 
UV photon (912 < X < 1 120 A) emitted by a star. The short 
wavelength cutoff is imposed by the usual assumption that 
enough neutral hydrogen is present to shield H2 from exci- 
tation by stellar Lyman continuum (Lyc) photons. The rate 
of molecular excitation is proportional to the stellar flux ab- 
sorbed by the ro-vibrational population of the ground elec- 
tronic state X L+. Fluorescence follows excitation, produc- 
ing a highly structured emission in the 912 - 1650 A band- 
pass, as electrons fall back into the excited vibrational lev- 
els of X '£+. Approximately 11 - 15% of the time, the flu- 
orescent pumping process leaves the molecule in the vibra- 
tional continuum of the ground state (X ), located 4.48 eV 

above the (7",v") = (0,0) level, 9 r esulting in its spontaneous 
dissociation into H(ls) +H(ls) dStecher & Williams! [19671 
iDalgarno & Stephensl [19701 IDraine & Bertoldilll996l) . Elec- 
trons ending up below the vibrational continuum cascade to- 
ward the ground vibrational state v" = via slow quadrupole 
transitions with radiativ e lifetimes ~ few days and longer 
(IWolniewicz et al.ll 19981). producin g the infrared fluorescence 
spectrum (iBlack & Dalgarnoi r976). 

Lyman and Werner band absorption of the central star con- 
tinuum unavoidably results in ultraviolet fluorescence, so why 
was it not detected? To answer this question we estimate the 
total brightness by simply calculating how many stellar pho- 
tons are absorbed by H2 at a given distance, which we take to 
be the transition zone, and re-emitted into 4tc steradians. This 
will give an upper limit to the total brightness of the H2 in the 
diffuse medium. We will neglect absorption by dust. 

Using our hydrogen absorption model we find 30% of the 
photon flux in the 912-11 20 A bandpass is absorbed. Reradi- 
ation into 471 steradians at 0.34 pc yields a total brightness of 
3.4 x 10~ 5 ergs cm~ 2 s _1 sr _1 (assuming an average energy 
per photon of 1.6 x 10~ n ergs). Neglecting the fluorescence 
redistribution into specific lines, and averaging over a 910 - 
1660 A bandpass, we find an upper limit to the mean bright- 
ness of the continuum fluorescence of 1 . 1 x 10 -18 ergs cm" 2 
s _1 A -1 arcsec~ 2 . This is a factor of ~ 50 below the back- 
ground estimate shown in Figure|3] So we find the continuum 
fluorescence is not very bright and shows no sign of contribut- 
ing to the ground state ro-vibrational population, as there is 
little in the way of absorption arising from vibrational levels 

9 Direct radiative dissociation of the X 1 state is dipole forbidden. 



v" > 2. This is in contrast to the strong indication of contin- 
uu m fluorescence fou nd in the reflection nebula NGC 2023 
bv lMever et al.l (1200 ll) . where they observed absorption lines 
from all v" < 14 states of X in a STIS spectrum of the 
central star HD 37903. 

The low levels of continuum fluorescence found here are 
consisten t with the infrared spectr oscopic diagnostics pre- 
sented by Zuckerman & Gatlevl d 19881) . who found the ratio of 
S(l) (1-0) / S(l) (2-1) (at 2.122 and 2.248 fjm respectively) to 
be ks 10, placing it in the thermal regime well above the va lue 
of 2 expected from fluorescence dBlack & Dalgarnol 19761) . If 
the S(l) (1-0) line were produced by fluorescence, then the 
10~ 4 e rgs cm~ 2 s _1 sr _1 measured by Zuckerma n & Gatleyl 
( 1988) would imply a total ultraviolet brightness ~ 850 times 
this value (in energy units, 50 times in photon units) and 
would have easily been detected by the rocket experiment. 

The detection of Lyoc fluorescence pumped through two res- 
onant transitions in the v"= 2 vibrational level of the X 

electronic state (R(6) (1-2) at 1215.730 A and P (5)(l-2) at 
1216.073) reported recently by iLupu et al.1 (120061) . serves to 
emphasize that in a radiation bounded nebula two-thirds of 
the Lyc photons em itted by the star are converted to Lyoc pho- 
tons (ISpitzerll998l) . The contrast between the Lya and contin- 
uum fluorescence is enhanced in our case, because the num- 
ber of Lyc photons exceeds the number of photons in the 912 
- 1120 A bandpass b y a factor of 15. 10 Even so, modeling 
by ILupu et al.l y006) shows the Lya fluorescence rate is too 
low to cause detectable deviations in a thermalized ground- 
state population at a temperature of 2000 K. We conclude that 
the observed ro-vibration population distribution of the hot H2 
is dominated by collisional processes. The dev iatio ns from a 
pure thermal distribution will be addressed in § 15.41 

5.3.2. Molecular Hydrogen Formation and Destruction 

The observed neutral atomic to molecular hydrogen ratio 
( j^pjg = 127) can be used to constrain formation and destruc- 
tion processes in the diffuse nebular medium using the for- 
mula, 

dn2 
dt 



= F-Dn 2 . 



(2) 



Ms the sum over all formation processes, D is the sum over 
all destruction processes, and «2 is the density of H2. We 
will not attempt to provide a comprehensive listing of all the 
processes, rather we will make some limiting assumptions to 
develop the constraints. 

In the limit where formation is absent (F = 0) the solution 
for ni(f) is an exponential with an e-fold time equal to D -1 . 
If we assume the destruction process is radiative dissociation, 
then D — XPdisJ^pump, where p^ s is the mean dissociation 
fraction, ^ pump is the mean probability per second of an up- 
ward radiative transiti on and X, the stand a rd int erstellar radi- 
ation field multiplier. IDraine & Bertoldj| (1 19961) have conve- 
niently tabulated K, pwnp and pdi SS for a number of transitions 
for % = L neglecting self-shielding. They find fairly narrow 
ranges for ^ pump and p c n ss and weighted mean values of these 
products over the H2 ground state ro-vibrational distribution 

10 Additional enhancement of Lya fluorescence results from the concen- 
tration of excitations into a single v'= 1 level of the B 1 state, as opposed to 
the continuum case where excitations to a multitude of vibrational levels, < 
v' < 20 for B 'Z+and < v' < 6 for C n„, compete for a relatively smaller 
number of photons. 
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are relatively insensitive to the actual distribution for our pur- 
poses here. We will use pdi SS = 0.15 and L, pw np = 3.5 x 10~ 10 
s _1 , yielding D\ = 5.3 x 10~ n s _1 as a representative disso- 
ciation rate for % = 1 . 

To compute x we use the Rauch model, which has a total 
photon luminosity of 1.2 x 10 45 ph s" 1 in the 912 - 1120 A 
bandpass at the stellar surface. At 0.34 pc, the photon flux is 
8.9 x 10 7 ph cm~ 2 s _1 , yielding a rather modest scale factor 
of % = 7.4 relative to the average ISM background of 1.2 x 
10 7 ph cm" 2 s- 1 dDraine & Bertoldilll996l) . The resulting e- 
fold time for radiative dissociation is 21 years and it would 



to lead to dissociative excitation. For that a radiative pro- 
cess is required such that H2{X l Y.+) + y ultimately results 
in H(l s) +H(3p), which can the n branch into Lyp or Ha 
+ Lya. IStecher & Williams! d 19671) discussed two possible ra- 
diative paths that lead to the dissociation of H2. In § |5.3.1| we 
described the familiar spontaneous dissociation process that 
leads to H(ls) + H(ls) following excitation into the Lyman 
(B £+) and Werner (C l Yl LI ) bands. The other path (with mul- 
tiple channels) is, 



only take 100 years for the ratio to reach its currently Hl ( X S s { 



N(H 2 

observed level. We conclude the source term for the diffuse 
medium F ^ 0, i.e. the H2 in the diffuse medium must be 
constantly replenished. 

We now turn to the equilibrium limit, i.e. ^ = 0. The for- 
mation term is F = yn«i where the neutral hydrogen density 
is n\. We consider either formation of H2 by dust, or by as- 
sociative detachment H +H~ — > H2 + e ~. In the for mer case 
Yd = 3 x 10~ 17 cm 3 s- 1 ( Natta & HollenbachlfT998h and n is 
the gas density. In the latte r case n is the electron den sity n e 
and y_ is given by Eq. 15 of iNatta & Hollenbachl d!998l) . The 
rate coefficient (y_) depends explictly on hydrogen ionization 
fraction and temperature. At T = 2000 K we find y_ rj y ( / 
when the ionization fraction x + = n + /{n\ +«+) = 0.5. When 
x+ = 0.05, y_ » 15y (/ . 

Equating formation and destruction we solve for the total 
density n required for formation to keep pace with dissocia- 
tion. With % = 7.4 the dissociation rate is D = 3.9 x 10~ 10 s -1 

and we find, n = Ul = 10 5 cm -3 for the dust and the x + = 0.5 

Y n\ 

7000 cm -3 . These den- 



H(\s)+H(nl) 



Pre — and direct 
Phot oionizat ion, 



- dissociation, 



(3) 



cases. Forx + = 0.05 the den sity is 
sities are consistent with the lMeaburn & Lopez! (1 19931) clump 
estimates for M27, but they are much higher than the gas and 
electron densities we found from the C I analysis in § 14.4.2 



for the transition zone or the electron density given by Barker 
( 1984) for the ionized region. 

We conclude that the density in the diffuse medium is not 
high enough to support formation of H2 either on dust grains 
or through the associative detachment processes and that the 
most likely source of excited Hg is from photo-evaporatin g 
clumps (c.f. iHuggins et all 120021: iLopez-Martfn et all 1200 lb . 
Whether the clumps are reservoirs of H2 that formed long ago 
or are currently active in forming H2 by these processes is 
an open question. The transport process that moves molec- 
ular hydrogen out of the clumpy medium and constantly re- 
plenishes the diffuse medium is a dynamical process, which 
results in the exposure of H2 to Lyc radiation. 

5.4. Beyond Spontaneous Dissociation - Lyc Dissociation 
and Photoionization 0///2 

We have two outstanding puzzles. One is the apparent "two 
temperature" slope in the ro-vibration population (Figure [9). 
The other is the excess Ha/Hp ratio in the absence of any sig- 
nificant stellar or nebular reddening. We speculate here that 
a confluence of environmental factors - molecular hydrogen 
interacting and excited by a ionized medium that is optically 
thin to Lyc radiation - conspire to contribute to the production 
of Ha emission at a few percent of the radiative recombina- 
tion rate and to create the two temperature slope. 

We know that while the energies of thermal electrons in 
the nebula are high enough to populate the X x £+ vibra- 
tional levels above v" > 0, they are not nearly high enough 



with n >2 and 1 = s, p, d ... the angular momentum of the ex- 
cited atom. The photoionization channel can further undergo 
dissociative recombination with low energy electrons, 



H+(X 2 Y+)+e- -^H(ls) +H(nl). 



(4) 



These channels are not usually considered to be important 
for H I regions (PDRs) because the photoionization, pre- 
dissociation and direct-dissociation processes all have thresh- 
olds above the Lyman limit. In such an instance the H2 is 
usually cold and opti cally thick neu tral hydrogen shields the 
H2 from dissociation (Spitzer 1948). However for H2 in H II 
regions it is not at all clear that this assumption is justified. 
At the interface of a mostly molecular clump embedded in an 
ionized medium, as we encounter in M27, there will be a zone 
where the H I and H2 are optically thin to ionizing radiation 
and the H2 is excited (hot). In this case the ionization and 
dissociation of H2 by photons with energies above the Lyman 
limit at 13.598 eV should also be considered. 

Direct-dissociation paths proceed through excited ro- 
vibrational levels in the B '£+, B' '£+ and C H u upper elec- 
tronic states of H 2 , while the B'"£+, D l U lt , D n U u and 
D" l Yl u states can under go pre-dissociation as well as direct- 
dissociation. The experimental distinction is that a pre- 
dissociation process shows structure (in the absorption cross- 
sections as a function of wave length) while direct-dissociation 
does not dLee & Judgdll976l) . The structure is a consequence 
of the perturbative coupling with the dissociative continuum 
of the overlapping excited electronic states. 

The vibrational continua of the upper electronic states con- 
verge into Rydberg orbitals, which are essentially n > 2 prin- 
cipal quantum number electrons in orbital around the cation 
(Hi" ) core. At asymptotic nuclear separation these orbitals be- 
come free //(Is) +H(nl). Transitions between X and its 
own vibrational continuum (H(ls) +H(ls)) at 4.48 eV can 
proceed only by quadrupole radiation, while transitions be- 
tween X 1 £+ and the upper singlet electronic states are dipole 

allowed. The vibrational continua of B '£+, B' '£+, and C l Yl u 
converge to H(ls)+H(2l) at 14.67 eV (844.8 A), D 1 n„ and 
B" *£+ converge to H(ls)+H(3l) at 16.56 eV (748.46 A), 
D n Yl u converges to H(ls) +H(4l) at 17.22 eV (719.74 A) 
and D" l U u converges to H (Is) +H(5l) at 17.52 eV (707. 175 
A), measured with respect to (/", v") = (0,0) of H 2 X The 
onset of photoionization (H2 + y — > Hj + e~) is at 15.43 eV 

(803.7 A) in between the H(\s)+H(2l) and H(\s)+H(3l). 
The right side of Figure [18] shows the energy levels of Ryd- 
b erg states in relation to the ca tion ground state. 

Glass-Mauiean et al. (1986) give cross sections for the 
photo-excitation of Lya (//(Is) +H(2l)) from the direct- and 
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FIG. 18. — Energy level diagram illustrating the threshold energy ro- vibration levels (J",v") that can be excited to Rydberg orbitals (H(ls) +H(nl),2 < n < 
°°,/ = s,p,d...) with photons energies < 13.6 eV. The threshold for excitation into the cation ground state is also shown. Normalized Maxwell-Boltzmann 
energy distributions for T = 10000, 4000 and 2000 K show schematically the relative change in the energy level population as a function of temperature. 



pre-dissociation of H2(X Zj ). The total cross section peaks 

at the threshold (14.67 eV), reaching 5 x IO7* 8 cm 2 and 
falls monotonically thereafter. Cook & Met zgerl d 19641) have 
shown the detailed structure in the H2 photo-dissociation and 
ionization cross section from 550 - 1000 A. The photoioniza- 
tion cross section peaks near 18 eV at w 10~ 17 cm 2 . Some 
of the discrete dissociation peaks to the red are « a fac- 
tor of 3 higher at a resolution limit of 0.5 A. iFlannerv et all 
d 19771) have calculated and tabulated cross sections for excita- 
tion out of vibrationally excited states of H2(X v") + y — > 

H^ (X 2 £+, v') + e~ as a function of wavelength, assuming AJ 
= 0. 



lYan et all d 19981) produced a piecewise continuous function 
for the photoionization cross section valid over the interval 
300 > Epf, > 15.4 eV. At high frequencies the cross section 
for H2 ionization is 2.8 times that for H I. Above 18.09 eV H2 
can undergo dissociative ionization, H2 + y — ► H + H + + e~, 
or even H + + H + + 2e~, but the branching ratios into these 
channels are sm all over the wav elengths of Lyc emitted by the 
central star (see lYan et alj 19981 and references therein) where 
th e yield for dissociative reco mbination should be p c i r w 1 . 

IStecher & Williamsl d 1967b have noted how the thresholds 
for transition into the upper level continua a re lower for ex- 
cited r o- vibrational (/", v") levels of X ICook & Metzgeil 
( 1964) observed this effect in the redward broadening of the 
dissociation thresholds by « 14 A in their experiments con- 
ducted at room temperature. It is actually possible to disso- 
ciate and photoionize H2 with photon energies less than the 
Lyman limit, for a high enough ro-vibrational level (J", v"). 
We quantify this in Figure [18] where the energies of the ro- 
vibrational levels in the X state in relation to the asymp- 
totically free states are shown. Horizontal lines drawn through 



the (v", J") energy levels indicate the thresholds for excitation 
with photons < 13.6 eV in energy. Maxwell-Boltzmann En- 
ergy distributions with temperatures of 2000, 4000 and 10000 
K are drawn to indicate schematically how the relative popu- 
lations of the upper levels change as a function of T. 

Rotational states above J" = 13 in v" = 0, J" = 9 in v" = 
1, and J" = 3 in v" = 2 can be excited to H(ls)+H(2p) by 
photons with energies < 13.6 eV. The steepening of the slope 
of the level population in Figure [9] might result from a radia- 
tive depopulation of the higher J" states into the asymptotic 
H(ls) + H(2p) state. The populations of the upper levels be- 
come selectively depleted, a process we refer to as "radiative 
skimming." However, the slope change in the v" = 0, 1 curves 
are observed to occur at at J" w 6, 8 respectively, which are 
smaller than given above. We take this as evidence for skim- 
ming by Lyc photons. 

We also see from Figure fT~8l that the creation of H(ls) + 
H(3p) without the benefit of Lyc photons would require ex- 
tremely high temperatures > 4000 K. The two step process 



H 2 (X 1 E+) + Y-.H+(X z j:+) + e 



H(ls) - 



ergetically more accessible than the direct H2( X £j ) + Y 



H(3p) is en- 

H(ls)+H(3p). Observations by |Lee_&Judgi il976|) of the 
direct production of B aimer series fluorescence from H2 by 
direct- and pre-dissociation support this conjecture. They find 
the dissociation yield with respect to total dissociation + pho- 
toionization cross section to be on the order of a few percent. 
For example, at 7 16 A the cross section for production of Ha 
was 0.23 x 10~ 18 cm 2 , which is about a factor of 40 lower 
than for the production of H^ . 

Hj is also homonuclear dLeach & Moss! IT995). so radia- 
tive transitions among the 18 vibrational levels of the X 2 £+ 
ground state are not dipole allowed. The vibrational contin- 
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uum of X 2 £+ (H + +H(ls)) is only 2.646 eV above the 
(/", v") = (0,0) ground state, so electron energy need not be 
high to produce dissociative recombination into the overlap- 
ping Rydberg states. Excitation to higher bound states, ( 2 £+, 

2 n„) have poor overlap integrals (Franck-Condon factors) 
yielding small branching ratios to these dipole allowed upper 
levels. Once created, H^ will undergo rapid dissociative re- 
combination with free electrons, H^ + e~ — > H(\s) + H(nl), 
producing a fluorescent cascade of Lyman and Balmer emis- 
sion. The amount of Ha and Lya produced will depend upon 
the vibrational distribution in the cation energy levels and the 
el ectron impact energies . 

Takagi (2002) has shown in detail how the dissociative re- 
combination cross sections of the individual H(nl) channels 
depends strongly on the vibrational level of the Hj X 2 £+(v") 
state. The cross sections are generally highest for e~ colli- 
sions with the lowest energies and fall in proportion to the 
inverse of the electron energy. Typical values near zero elec- 
tron energy range from 10~ 15 - 10~ 17 cm 2 ignoring rotation. 
Including the rotational state introduces resonances and can 
strongly modulate the cross sections at low electron energies. 
ISchneider et alj ( 1 19941) have calculated total dissociative re- 
combination rates for H^ for hydrogen plasmas with tempera- 
tures between 20 - 4000 K for individual ro-vibrational states 
of Hj , but they provide no estimate of the branching ratios 
into the individual H(ls) +H(nl) final states. 

5.5. Optically Thin Photoionization Destruction of Hi 

We have suggested that photoionization by Lyc photons of 
excited H2 in an optically thin medium ultimately results in a 
fluorescent cascade of Lyman and Balmer emission. Verifica- 
tion of this hypothesis will require exploring the detailed bal- 
ance of hydrogen ions, atoms and molecules in the recombi- 
nation regime between Case A (the optically thin limit where 
Lyman series emissions are not reabsorbed in the nebula) and 
Case B (the optically thick limit where Lyman series emis- 
sions are reabsorbed on the spot). 

In the previous section we have given references to the nec- 
essary molecular cross sections for computing the photoion- 
ization rates out of the vibrationally exc ited molecule and into 
the excited cation dFlannery et al.lfl977l) . along with the cross 
sections for electron impact diss ociation of e xcited cation 
states into the Rydberg manifold (Takagi 2002). The calcu- 
lation of the Lyman and Balmer series emissivity resulting 
from the two step process of molecular photoionization fol- 
lowed by dissociative recombination goes beyond the scope 
of this paper. 

Here we show an example of the competitiveness of pho- 
toionization with spontaneous dissociation for the destruction 
of H2 in the optically thin limit. The ionization rates at the 
transition zone (0.34 pc) for H2 and H I were calcul ated by in- 
tegrati ng the cross sections for th ese species, using lYan et alj 
(1998) for H 2 and ISpitzeri (fl998l Eq. 5-4) for H I, over the 
SEP of our stellar mode l from A, < 911.7 A according to 
Spitzer (Eq. 4-36, 1998). We do not include contributions 
from direct and pre-dissociation processes, which are likely 
to be substantial for excited H2, so our dissociation estimate 
will be conservative. We find C,H2^H2+ = 3.4 x 10~ 9 s _1 
compared to C,h^h+ = 1-4 X 10~ 9 s -1 . The ratio is 2.4, so 
we conclude that in the optically thin limit it is easier to ionize 
the hydrogen molecule than the atom. 

In an optically thin medium at the inferred radius of the 
transition zone the molecular photoionization rate t s H2^H2+ 



exceeds the far-UV pump rate X^pump (§ 15.3.21 by a factor of 
1.3. If we assume that the yield for dissociative recombination 
(Hj + e~ — > H(ls) +H(nl)) following creation of the cation 
is 100%, then the photoionization mediated dissociation of H2 
can be expected to be Pdr^H2-^m+ I (PdissX^pump) ~ 9 times 
the spontaneous dissociation. This process increases the rate 
required for the clumps to resupply the diffuse medium with 
H 2 in the steady state limit. We note that lNatta & Hollenbachl 
( 1998) do not include Lyc dissociation of H2 in their model. 
This may be why they find the shell mass is dominated by H2 

as opposed our finding that 3> 1 ■ 

6. CONCLUSIONS AND SUGGESTIONS FOR FUTURE 
INVESTIGATIONS 

M27 is an excellent laboratory for testing theories of for- 
mation, excitation and destruction of molecules and dust, be- 
cause the stellar temperature, mass, gravity, as well as the neb- 
ular distance, masses, abundances and excitation states, are 
well quantified. 

We find that the diffuse nebular medium is generally inhos- 
pitable to molecules and dust. The diffuse H2 accounts for a 
small fraction (~ 0.3 %) of the total nebular mass. H2 is easily 
destroyed by the stellar radiation field in the diffuse medium 
on a timescale ~ 20 years. A steady state abundance requires 
a source, most likely the clumpy medium. If the clump den- 
sity is high enough, they may support molecular hydrogen for- 
mation via dust or mediated by radiative detachment, although 
we cannot rule out photoevaporation of relic material. 

The finding of neutrals and molecules at a transition veloc- 
ity between slow-high and fast-low ionization material sug- 
gests collisional interaction with charged particles is the most 
important process for establishing the ob served ro- vibration 
level populations in molecular hydrogen. iLupu et all d2006l) 
finds the rate of Lya pumping is not high enough to change 
the level population of a thermal distribution and we find the 
rate of continuum pumped fluorescence is even lower. How- 
ever, if the neutral hydrogen is optically thin in the transition 
zone, then stellar Lyc can radiatively skim the upper rotational 
states of v" = 0, 1 and create the apparent break in the high and 
low temperature slopes at the observed levels, but this is likely 
to be a second order process. 

Dissociation and ionization of the hydrogen molecule by 
Lyc photons can also lead to the production of Lyman and 
Balmer series emission and may provide an explanation for 
the excess Ha/Hp ratio observed in the apparent absence of 
extinction by dust. In the optically thin limit Lyc mediated 
destruction rates can exceed by an order of magnitude the 
spontaneous dissociation rate and should be considered when 
calculating the equilibrium abundance of excited H2 in H II 
regions. 

At this late stage in the evolution of M27 most of the orig- 
inal AGB atmosphere has been dissociated and accelerated to 
velocities in excess of 33 km s but not above 68 km s . 
A high speed radiation driven stellar wind is not detected in 
the present epoch. The mechanism for acceleration remains 
unidentified. 

The analysis of the formation and destruction of H2 rests 
in part on the determination of the H I velocity structures 
identified in the Dwingeloo survey data to guide the column 
density determination carried out with the STIS and FUSE 
data. However, our conclusions regarding the overall abun- 
dance of molecular hydrogen with respect to H I and the need 
for a molecular source for the diffuse medium are immune to 
changes in the atomic abundance by factors of 10. Neverthe- 
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less it would be prudent to check the velocity structure found 
here with a high resolution mapping using a telescope with a 
smaller beam (« 4'). It appears that the velocity separation 
of M27 permits discrimination with respect to the galactic H I 
background. This and a reasonable signal strength have al- 
lowed the apparent detection o f H I in a PN where on ly upper 
limits have existed previously (Sc hneider et al.|[T987l) . 

Our lowest bound for N(CO)/N{H 2 ) < 1(T 3 in the dif- 
fus e nebular medium, is a factor of 30 higher than ratio used 
by Huggin s et al.l 0996) to constrain the CO clump mass. 
Should a HST servicing mission 4 be successful in resurrect- 
ing STIS then higher signal-to-noise data could be acquired to 
examine the N(CO) /N(H2) more carefully, using a number of 
rovibrational bands. It should be possi ble to place limits on 
logJV(CO) < 13 (e.g. lBurghet"aill2007h . 

We emphasize that our analysis applies to the direct line 
of sight gas, the diffuse nebular medium, which we find to 
be dust free. The existence of clumps with column densi- 
ties of CO in excess of 10 16 cm~ 2 immediately adjacent to 
regions with CO column densities < 10 14 cm~ 2 attests to 
the steep density gradients within the nebular environment 
and the probable low covering fraction of the densest clumpy 
medium. The clumps may well be reservoirs of material de- 
pleted from the diffuse medium. Measuring the abundance 
of metals on either side of the neutral transition velocity 
could provide information on whether photo-evaporation of 
the globules is an important pr ocess for the enrichment o f 
metals in the high velocity zone dMcCandliss & Krukll2007h . 
In addition, infrared spectral and imaging observations with 
Spitzer could tell us much more about the prevalence of dust 
and molecules elsewhere in the nebula. Mapping of the neb- 
ular Balmer emission line ratios with respect to the H2 S(l) 
(2-1) / S(l) (1-0) emission ratio would be especially useful 
for further exploring the connections between H2 production 
of Balmer emission and the prevelence of thermal or fluores- 
cence processes in the clumpy medium. 

Dur ing the review of this paper, work by lO'Dell et al.l 
(l2007h was brought to our attention. O'Dell et al. derive 



the total flux in the infrared H2 emission lines in the Helix 
Nebula and argue that only Lyc radiation from the star has 
enough power to sustain the H2 IR emission at the observed 
levels. They also propose that charge exchange between H2 
and + could be a significant dissociation channel. Future 
work examining the efficiency of their proposed charge ex- 
change mechanism for mediating the dissociation of H2 in 
comparison to the Lyc dissociation process proposed here, 
would be a useful addition to the nascent field of Lyc dom- 
inated PDR's. 
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